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*  mtogccnox 

lafcrwaoQ  abcvii  the  $urti-’o  of  th*  ci:th  has  for  long  been  ootained  by  a  vide 
variety  of  methods.  With  the  ability  to  ascend  above  the  sort  are  baa  occurred  what  is 
now  referred  to  as  'reaote  sensing’  tsee  Kef  1),  the  acquisition  cf  inf onset too  frees  a 
distance ,  first  by  aerial  phetegruphy  and  latterly  by  ea"-th  satellite  observation. 

Observation  of  the  earth  from  satellites  has  been  aide  by  photography  and  by  other 
imaging  net  hods.  In  soa  cases  the  photographs  can  c«  physically  returned  to  earth 
(cr  Sky  lab),  but  this  is  cot  always  cuavenient,  and  is  certainly  not  so  for  satellites 
designed  to  make  cony  observations  over  a  long  period  of  tine.  As  an  alternative  to 
cooventioaal  photography,  therefore,  the  earth's  surface  nay  be  optically  scanned  by 
sane  nesns  and  the  observations  converted  to  a  fora  which  can  be  transmitted  to  earth. 

In  this  case  the  data  can  be  reconstituted  into  a  tors  which  resent  I es  a  photograph,  mad 
which  will  be  n 'erred  to  in  this  Report  as  a  'picture*.  This  Report  is  concerned  with 
digital  data  and  pictures  provided  by  the  mull i-spectral  scanner  (ICS)  on  the  Landsat 
series  of  earth  resources  satellites. 

Images  provided  by  Landsat  S5S  cover  an  area  approximately  I  S3  kn  square,  viewed 
vertically  downwards,  and  are  of  tec  immediately  recognisable  by  a  viewer  whe  is  familiar 
with  the  portion  of  country  depicted.,  Recognition  is  usually  by  means  of  salient 
features,  such  as  the  shape  of  coastlines  or  the  position  of  roads.  However,  a  small 
amount  of  spatial  distortion  is  present,  sc  it  the  precise  location  of  any  particular 
feature,  or  any  exact  distance,  is  required,  it  is  usually  necessary  to  resort  to  inter¬ 
pret  iw  techniques.. 

Perhaps  the  mest  convenient  treatment  ter  Landsat  MSS  data  is  to  convert  or 
‘transfer**  it  to  a  known  cartographic  coordinate  system.  The  information  is  then 
presetted  in  a  fora  which  is  familiar,  whose  dimensions  (scale  and  other  properties)  are 
understood,  which  nay  include  calibration  grid  lines  and  which  nay  be  directly  cmparid 
with  a  conventional  map  of  the  same  scale,  to  establish  locations.  Satellite  imagery  may 
then  be  studied  by  a  wide  variety  of  users,  without  their  having  to  be  concerned  with  the 
nature  and  details  of  the  recognition  and  transformation  process. 

Two  distinct  forms  of  iaroraatio  are  in  general  available  for  precise  dimensional 
evaluation  of  photographs  or  photogra;  a-like  images.  One  form  of  information  relates  to 
the  details  of  the  imaging  systen,  for  example  in  microphotography  the  magnification  of 
the  system,  and  in  telephotography  knowledge  of  any  distortions  introduced  by  the  optical 
system  as  instanced  by  the  difference  between  the  behaviour  of  a  sliding-back  and  a 
conventional  camera.  The  other  fora  of  information  is  knowledge  about  the  object  being 
photographed,  for  example  a  measuring  rod  may  l.  included  in  a  photograph  of  an  archaeo¬ 
logical  dig  to  provide  the  local  sea’e. 

The  actual  process  of  Landsat  image  identification  and  transformation  makes  use  of 
both,  types  of  knowledge,  c  those  concerning  the  imaging  system  and  concerning  the 
earth's  surface.  Certain  features  of  the  image  can  be  'recognised*  and  these,  together 
with  a  knowledge  of  the  geometry  of  the  optical  imaging  process,  allow  the  remainder  of 

the  image  to  be  fitted  ir.to  its  correct  place. 


The  recognisable  terrestrial  features  are  referred  to  as  'ground  control  points* 
igcpi) ,  and  to  he  usable  they  must  he  identifiable  both  on  the  satellite  inage  and  also 
on  soar  suitable  calibrated  map. 

I,!  The  Landsat  MSS  image 

The  Landsat  KSS  is  described  in  detail  in  Ref  Z.  Briefly,  the  satellite  acrm  in 
a  near  polar  orbit  (inclined  at  about  BO  degrees  to  the  equator)  at  a  nominal ly  constant 
height  (about  ?OC  ta)  above  the  earth's  surface.  As  the  satellite  travels,  the  surface 
vertically  below  it  is  scanned  in  a  cross-track  direction,  so  that  the  observations  form 
a  raster  scan  of  the  earth's  surface.  The  inage  is  resolved  into  picture  eleaents  or 
pixels  of  about  80  a  square,  witn  soae  overlap  of  pixels  along  the  scan  line,  so  that  each 
pixel  is  about  b0  a  ir  that  direction.  Observations  are  made  at  four  spectral  wave  lengths, 
but  this  need  not  be  considered  here  except  to  note  that  in  practice  the  nost  convenient 
waveband  for  determining  gcps  is  the  so-called  band  7,  (0.8  to  1.1  microns  wavelength, 
in  the  near  iafra-red. 

From  this  description,  sane  properties  of  the  resultant  image  may  be  deduced. 

Firstly,  the  scale  of  the  image  in  the  direction  of  satellite  travel  should  be  constant, 
since  tho  near  circular  o-bit  provides  a  near  constant  speed.  Secondly,  the  scale  in  the 
cross- track  direction  depends  upon  the  nature  of  the  scanning  mechanism  and  the  curvature 
of  the  earth's  surf  are.  (Variations  in  height  of  terrestrial  features  nay  be  ignored 
since  their  effect,  particularly  for  the  low  altitudes  of  the  majority  of  Britain,  is 
marginal  if  not  undetectable.)  Since  the  natures  of  these  twc  scan  effects  are  fairly 
accurately  krovr.,  tifjr  can  he  ccopeosated  for.  Giber  uawit  effect*  such  as  the  rotation 
of  the  earth  beneath  the  satellite  can  also  be  allowed  for.  The  resulting  image, 
referred  to  as  'system  corrected*,  is  in  the  form  of  an  oblique  Mercator  projection: 
indeed  Colvocoresses^  has  proposed  the  use  of  'Space  Oblique  Mercator*  projection  for 
Landsat  imagery. 


For  many  purposes  it  would  be  desirable  to  transform  the  image  to  a  more  conven¬ 
tional  map  projection,  « j  soae  form  of  Transverse  Mercator,  and  this  cm  in  principle 
be  done  by  computer  provided  the  latitude  and  longitude  of  some  locations  of  the  image 
are  accurately  known. 


In  practice,  certain  difficulties  arise,  and  these  stem  from  the  lack  of  perfection 
of  the  satellite's  movement .  The  orbit  is  much  that  the  height  shove  the  earth's  surface 
varies  by  approximately  tit  and  the  vehicle  itself  is  subject  to  muvemnts  in  all  three 
ot  its  rotational  axes,  ie  it  pitches,  rolls  and  ymrs.  All  of  these  movements  are  to 
seme  extent  measured  and  known,  but  net  with  sufficient  accuracy  to  allow  the  imp  to 
be  accurately  calibrated.  A  reasonable  level  of  accuracy  might  be  defined  Such  that  mo 
pixel  is  (ut  of  place,  with  reference  to  the  map  system,  by  more  than  its  own  dimensions. 

1.2  Geometric  trmtsfomation  of  Landsat  i— -ft-* 

Since  it  is  not  possible  to  relate  an  image  to  the  earth's  surface  with  sufficient 
precision  from  available  knowledge  of  the  operation  and  movements  of  the  1  ■«»««>  imaging 
system,  resort  has  therefore  to  he  made  to  knowledge  of  the  aubject  being  imaged.  This 


o?r 


■ty  b«  done  with  the  aid  of  geps  whose  posit  iot.  can  be  ikurmiud,  both  an  the  iss|f  fltd 
also  oo  the  earth's  surface.  Since  terrestrial  features  have  been  surveyed  in  great 
detail  in  Britain,  it  is  ceavericat  to  use  naps  for  the  accurate  determination  of  the 
position  of  geps  on  the  earth's  surface. 

When  sufficient  geps  have  beer,  located,  both  on  the  image  and  oo  a  nap.  acme  form 
of  mathematical  transformation  can  be  used  to  enable  all  points  in  the  image  to  be  trans¬ 
formed  to  the  map  coordinate  system,  this  process  being  referred  to  as  geometric  trans¬ 
formation. 

The  identification  of  geps,  both  on  an  image  and  on  a  map.  is  is  the  first  instance 
an  operation  which  is  done  by  human  beings,  as  no  suitable  pattern  recognition  system 
has  yet  bees  set  up  to  cocpar<  ar.  image  with  a  map .  (It  would  sees  feasible  to  perform 
this  operation  automatically  by  cogtiur,  for  certain  types  of  geps,  for  example. 

Selected  land-water  boundaries. 1  However  the  Lands  at  system  is  a  continuing  operation, 
so  repeated  images  are  produced,  and  it  has  proved  feasible  to  use  a  computer  to 
recognise  and  relocate  known  geps  in  these  repeat  images. 

A  distinction  should  therefore  be  clearly  drawn,  between  (a)  the  first  identifi¬ 
cation  and  location  of  geps  both  oa  an  image  and  on  a  map,  and  (b)  the  relocation  of 
known  geps  in  other  inages.  (Relocation  on  a  map  is,  of  course,  not  needed.)  This 
Report  is  mainly  concerned  with  (b) ,  and  describes  an  almost  entirely  automatic  computer 
process,  the  program  CCF.FIXD,  lor  relocation  of  specified  geps  within  a  repeat  image. 

The  subject  of  geometric  transt craat ioa  has  been  discussed  in  the  literature,  one 
cf  the  more  detailed  expos! lions  being  that  of  Shlim'*. 

I  .  3  Some  note*  oa  the  program  GCP.FrXD 

Jtany  countries  have  their  own  satellite  data  pro-.essing  facilities,  so  this  work 
has  been  primarily  concerned  with  images  of  Britain,  sc  the  mainland  and  surrounding 
islands,  but  excluding  Northern  Ireland  and  Ireland.  However  the  possibility  of  use 
elsewhere  has  been  considered,  and  those  features  of  the  work  relevant  only  to  Britain 
have  been  noted.  It  should  therefore  be  possible  to  adapt  the  program  QCP.FIND  for 
operation  in  other  countries,  with  fairly  limited  changes. 

The  problem  o-  relocating  geps  ir.  images  is  essentially  a  practical  one,  and  has 
been  treated  as  such  in  this  work.  The  processes  which  are  described  is  this  Keport  are 
ones  which  have  been  found  by  experience  to  work  satisfactorily  for  an  adequately  large 
proportion  of  geps,  on  a  vide  variety  of  images  of  Britain,  the  whole  operation  requiring 
a  reasonably  small  asvxiut  of  computer  facilities  and  operating  time  compared  with  other 
operations  performed  on  the  image.  The  methods  described  are  not  necessarily  the  'best* 
or  the  most  elegant  and.  indeed,  considering  the  variety  of  imiges  to  be  processed,  it  is 
doubtful  whether  a  'best*  method  exists. 

4  FIBST  LOCATION  OF  GfcXXB  CONIKH.  POINTS 


Since  the  first  identification  and  Ip  tioci  of  geps  aust  be  dome  before  their 
relocation,  the  former  process  will  be  des.  'bed  first.  The  author  of  this  Report  has 
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not  perforate  J4Y  of  this  Ic^ttign  work,  direct ly  niastit,  but  has  Hdt  several  suggestion 
which  have  proved  of  assistance  in  doing  it,  and  these  will  be  atstiootd  in  this  section. 

2. 1  Ground  control  point  library  systeg 

The  f  irst  gep  work  done  in  Space  Department  SAE  is  that  described  by  Williams^  who 
selected  about  >-  locations  ia  Landsat  inage  path  219  row  24  of  2  July  1977.  Subsequently, 
others  at  U£  have  selected  sufficient  geps  to  cover  the  whole  of  Britain.  Ceps  have 
been  identified  on  Lands it  inages  and  on  Crdrance  Survey  I : 500 X)  saps,  locations  being 
referred  to  the  National  Crid.  Since  a  correspondence  between  inage  and  nap  is  required, 
the  procedure  results  in  the  p-ovisioc.  of  two  pairs  of  values  for  each  gep:  the  column 
and  row  siab«r  of  a  pise*  in  the  relevant  Land  sat  image,  and  the  corresponding  National 
Crid  eastings  ar.J  northings,  recorded  to  an  aevuray  cf  about  0.2  tm  on  the  map,  for  the 
centre  of  the  selected  pixel,  oetweec.  SC  and  IOC  geps  are  selected  for  each  image,  -end 
for  the  whole  of  Britain  this  has  involved  the  selection  of  about  1200  geps. 

The  main  object  of  this  work  is  to  create  a  number  of  data  sets.  One  of  these  is 
a  ’master  sup-location  file',  tv  a  list  of  ail  geps  together  with  their  Bational  Crid 
locations.  Other  data  sets  are  created  for  each  image,  containing  a  list  of  the 
50-102  geps  tor  that  image  with  their  col warn  and  row  ntabers  in  that  image.  These  data 
pets  are  subsequently  used  to  create  a  set  of  gep  'chips*  as.  'ascribed  in  section  3.1. 

It  may  be  acted  that  for  Britain  it  is  only  necessary  to  process  alternate 
satellite  tracks,  as  there  is  a  full  SCI  overlap  of  adjacent  images.  Consequently  all 
of  Britain  can  be  examined  in  about  18  images. 

Sow  that  the  whole  of  Britain  has  had  geps  selected,  the  process  described  as  the 
main  subject  of  this  Report  allow.*.  all  other  images  to  have  their  geps  located  almost 
entirely  automatical ly .  Tfcus  the  ;imy  hundred  image:  with  usable  portions  of  Britain 
can  be  processed  rapidly  and  alm?s  automatical ly  as  required. 

2.2  Ground  control  point  local  ioc  i-o  maps 

Originally,  gep  locative  was  done  by  Villiam*  as  follows.  Ceps  were  visually 
identified  or.  image  and  map,  for  exaopl*  road  intersections,  recognisable  curves  of 
rivers,  shapes  cf  lakes.  For  each  such  ’eater*  the  nearest  image  column  and  row  were 
recorded,  together  with  National  Crid  easting  and  northing.  Subsequent  re— exmuuation  of 
these  geps  has  shows  that  the  locations  weie  not  always  as  accurate  as  might  be  wished, 
for  it  mgile  since  a  relative  rotation  of  about  fifteen  degrees  ex-sts  between  image  and 
map.  the  exact  location  on  the  curve  of  a  river  could  be  appreciably  in  error.  Again,  a 
road  intersection  could  be  nl* taken  for  an  adjacent  one. 

Following  experience  with  correlation  of  patches  of  images,  the  present  author 
suggested  that  small  patches  of  images  be  prepared,  to  the  asms  scale  as  the  maps,  and 
on  a  craospareat  base,  so  that  each  patch  could  be  laid  over  the  map  and  moved  about 
until  the  best  tic  was  obtained.  Furtlxsr,  since  the  angle  between  image  and  map  is 
known  to  at  least  an  accuracy  of  one  degree,  (section  7.2)  the  transparency  could  include 
a  line  to  indicate  the  relevant  east-west  direction  of  the  map.  This  line  could  then  be 
held  parallel  to  the  map  east-west  direction  and  the  fit  of  image  to  sq  would  then 


ltivc  1  w  only  ttdnsUtwn  in  lw-  directions,  vt'.hv-jt  r>U:ioc,  Finally,  to  aid  in  the 
tie  terminal  ion  ot  the  centre  ot  the  centre  pixel  ot  the  p*‘.h,  a  small  circle  could  be 
incorporated  during  the  preparation  o:  the  transparency.  These  aids,  together  with 
others  determined  in  the  light  ot  experience,  have  enabled  geps  tc  be  located  extremely 
accurately,  so  much  so  that  deficiencies  within  the  Image  Itself  have  been  exposed. 

These  deficievcu  s  relate  to  the  several  forms  of  system  correction  applied  do ring  the 
preparation  of  the  computer  compatible  tapes,  and  are  described  in  Ret  t.  Add  it  tonally , 
the  occasional  absence  of  scar,  lines  from  ar.  image  car.  be  cetected.  For  greatest 
accuracy  it  may  be  aecessarv  to  start  with  ’r*'  imagery*. 

2.)  The  transformation  matrix 

When  a  sufficient  number  ot  g._p»  have  been  Uts'.i-i  or.  bo"h  image  and  map,  the 
information  may  be  used  to  calculate  a  * trar.sfoim.it  ton  setrix*.  vhi.h  enables  all  image 
locations  to  be  converted  to  map  reference*,  or  vice  ver>a.  This  has  beer,  described  by 
Williams^  and  by  Shliea'*.  Various  typer  of  matrix  tar.  be  made,  according  to  the  ember 
and  accuracy  of  geps  xvailable,  and  the  accuracy  of  th«.  required  t  rons  format  ion,  A 
first-order  matrix  allows  the  iuage  to  be  rotated,  iraslite:  an-  sreareu  enty,  whereas 
a  higher— order  matrix  enables  istre  complicated  '.ransfjrauior.s  tc  be  performed,  such  as 
correction  for  non-linearity  of  the  satellite  scan  mechanism.  Preteat  eig>er:eaee 
indicates  that  ^  first-order  matrix,  gives  useful  results  when  applied  to  system 
corrected  data,  but  tha'  ar.  appreciable  improvement  can  be  obtained  with  a  higher-order 
matrix.  Since  'he  latter  would  also  b»  capable  of  handling  the  scan  lice  non-linearity 
present  in  noa-system-corrected  (cr  'raw')  data,  it  seems  likely  that  this  would  form 
the  optimum  system  for  highest  a.  cur  icy .  Other  aspv  t*  or  the  of  raw  data  will  be 
touched  on  later  in  this  Report. 

2.1  Xuadmring  system  for  ground  control  points 

The  whole  of  Britain  Vas  defined  in  section  l.i)  is  napped  by,  amongst  others, 

204  sheets  of  the  Ordnance  Survey  l:id.vW  Second  Seiies.  ta.b.  gep  has  been  given  a 
code  number  consisting  of  two  integers,  the  first  being  the  number  or  the  OS  map  on 
which  it  occurs  and  the  second  act  arbitrary  serial  number  cr.  that  map.  Sometimes,  doe 
to  overlap  of  map  sheets,  a  gep  e-curs  or.  m>re  than  one  sheet,  and  in  such  cases  it  is 
customary  tc  use  the  lower  or  lowest  ot  the  sheet  makers.  Failure  to  observe  this  rule 
carries  no  adverse  consequences  in  the  operation  of  the  coaputer  program  but  may  cause 
inconvenience  when  referring  to  maps.  Provision  is  made  for  serial  numbers  from  I  to  99 
inclusive,  although  the  average  number  is  less  than  six.  This  numbering  system  is 
described  here  mainly  so  that  the  detailed  operation  of  the  computer  program  may  more 
easily  be  understood. 


3  RELOCATION  OF  A  GftXTtD  CONTROL  POINT 

The  operations  described  in  section  2,  the  establishment  of  a  first  set  of  gcp«, 
have  been  performed  manually,  and  the  results  have  been  entered  into  data  files  ready 
for  use  by  a  computer.  The  determination  of  gep  locations  in  subsequent  images  esc  be 
dome  almost  co^letely  automatically  vSifc  the  aid  of  the  computer  program  GCP.FIXD,  which 
i'i  the  main  subject  of  this  Report.  This  program  makes  use  of  a  method  of  relocating 


ind  iv  id-tal  5  ps,  this  being  don*  by  the  subrowtir.*  UX'GCF  vlocat*  gep),  which  is 
described  ir.  detail  iu  u-.>r.  *.  £1*  nlUnisj  stauT.s  descr.be  how  Jill  gcps  in  ar. 

IM««  arc  located. 

3. 1  Ground  control  point  ’chips’ 

the  ta>k  ot  f$:a:  1  is-mg  the  location  of  a  g.p  in  an  image  is  essentia'  ij  one  of 
pattern  recognition.  Any  gep  is  by  definition  tie  description  oi  some  terrestrial 
feat-re  of  ttot»-z*ro  sire.  For  exiapte,  the  definition  “the  intorsec  t  ion  of  the  M-*  and 
5t>  a  t.'tv**"  whilst  it  rerers,  as  a  location,  to  the  point  where  the  centre  lines  of 
those  toads  inters*-*,  ir.  fart  only  has  meaning  with  reference  to  the  existence  cf  the 
roads  themselves  for  setse  short  distance  at.  and  to  either  side  of,  the  intersection. 

The  actual  location  cf  a  £ ;p  tier*.:  oro  ,  is  ie terxi-e;  by  the  desrription  ofc.k  particular 
point  in  its  s— iroiridine  context.  Tte  g.p  can  og t y  be  locs'e-d  by  iocatir.^  a  surrounding 
area. 

For  the  first  location  the  pattern  recognition  problem  is  that  of  relating  a 
portion  c:  ar.  image  with  a  portion  of  a  nap.  In  the  case  of  1-aaisa  >SS,  the  inage  is 
in  the  form  of  a  raster  of  rectangular  pixels,  whi.st  the  maps  used  '.OS  l : 30CC0)  consist 
of  lines  and  coloured  areas  printed  or.  papet  ,  :c  in  ntxi-digital  fora-  It  is  possible 
to  devise  various  methods  of  rela'ittg  image  and  map,  and  it  sertxc  th»t  the  method 
described  ir.  section  2  gave  probably  the  bei  t  ccsfeina' i.*:  of  accuracy  with  spe*  i  for  the 
one  thousar-c  cr  so  g.ps  repaired.  As  noted,  this  method  used  a  compute  r  for  the 
preparation  of  the  photographic  patches  of  image  which  were  coapa red  with  the  maps.  It 
seems  liwoly  that  a  vr»  automatic  method  cl  relating  im»ge  to  map  woal-.  aeed  considerab^.- 
aore  effort  to  develop  :har  the  method  describei. 

For  the  problem  of  relocating  geps  the  situation,  is  greatly  eased,  as  it  reduce*  to 
the  matching  o:  portions  of  oc.%.  image  with  another  image,  that  is,  both  elements  are  ir. 
the  same  format,  which  is  already  a  digital  one. 

In  readiness  tor  the  aatvhing  process,  small  'chips'  of  the  original  image  are 
i 

prepared.  These  chips  are  subimages,  cop'ed  fro®  the  original  image,  and  consist  cf 
19  rows  of  W  col-car.*  ot  pixel  values,  centred  or.  the  gep  nominal  column  and  row 
location.  A  'library*  of  gep  chips  has  been  prepared  and  these  are  available  for  use  by 
the  program  CCF.FlXo.  Fach  chip  forms  its  own  small  computer  data  file  and  the  file 
name  is  based  on  the  gep  number,  described  in  section  2  .  Thus  gep  number  172,3  has  a 

corre spending  chip  with  tile  name  Cl  72.3  in  a  specified  'user  area’  of  the  disc  data 
storage.  Chips  were  prepared  in  batches,  one  batch  per  image,  using  the  purpose-written 
program  IX.CHIPIX. 

The  problem  of  relocating  geps  is  made  much  easier  than  the  general  pattern 
recogsitioc  prooles  when  it  is  noted  that  (a)  the  satellite  height  varies  very  little, 
hence  the  size  of  a  gep  is  almast  identical  in  every  image,  and  lb)  the  direction  of  the 
satellite  path  (ie  its  angle  to  the  local  meridian)  varies  very  little,  hence  the  gep 
is  always  orientated  at  about  the  same  angle.  This  the  pattern  recognition  problem 
be  cases  the  manageable  oru_  of  co^aring  and  relating  images  of  almost  identical  size  and 
orientation.  It  is  true  that  the  object  being  imaged  may  change,  from  one  satellite 
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pass  to  another,  a  field  say  be  ploughed  cn  one  occasion  and  have  growing  crops  on 
another,  and  such  changes,  discussed  later,  can  sake  the  relocation  problen  sore 
difficult,  if  not  impossible  in,  for  example,  , he  cose  of  a  stwc* cover* i  rcene,  or  of 
course  a  c loud-covered  gcp. 

In  practice,  it  has  been  found  that  gcp  relocation  <%n  usually  be  satisfactorily 
perforurd  by  pattern  matching  a  ‘chip*  fro*  one  image,  with  another  image,  as  described 
be  low. 

3.2  Pattern  matching 

Ihe  nethod  of  establishing  the  correspondence  of  chip  and  new  image  location  is  by 
means  of  the  correlation  coefficient,  “C  ,  calculated  between  the  !9  by  !9  b»  361)  pixels 
of  the  chip,  and  a  19  by  19  pixel  section  of  trie  new  image.  The  correlation  coefficient 
is  calculated  for  ary  image  location  by  the  convent  tonal  mfthod: 
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where  x. 


ith  pixel  value  of  the  chip 
y -  “  ith  corresponding  value  for  the  image 


»  average  of  all  361  chip  pixel  values 
y*  v  average  of  all  361  image  values 
and  the  summations  are  coalucted  tor  i  *  I  to  361  , 
Expression  U)  may  be  rearranged  in  the  fera 
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where  N  •  361  ,  < he  total  r.usber  oi  x-  or  y  values,  and  the  summation  is  as 

before. 

Expression  (2)  i*  more  convenient  tran  t I 2  for  computer  evaluation  because  it 
does  not  involve  the  pn— calculation  of  the  average  values  x^  and  y^  .  In  subroutine 
LOCUCF  r.he  calculatic>n  i;  seen  ir.  the  lorn 


to 


bincc  Rf  only  contain*  x^  values,  it  guy  bv  calculate  one*  only,  for  the 
particular  chip,  a  J  will  then  be  available  for  all  further  correlation  calculations  fo 
that  chip,  whereas  the  terns  R!  and  R)  m»st  be  recalculates  fer  each  position  of  tlie 
chip  in  relation  to  the  new  iuj<. 

I*  either  t  '=e  image  or  the  chip,  or  both,  should  be  uniform  over  their  19  by  19 
extent,  then  the  calculation  ct  R!  produces  the  value  rero.  In  the  case  of  a  uniform 
chip,  then  Rl  is  also  zero,  and  for  a  ur.iform  image  patch,  Rj  is  aero.  Ir,  either  case 
the  CC  calculation  is  ir.de  te  ruinate ,  be  ing  the  result  of  zero  divided  by  tero.  An 
attempt  to  divide  by  aero  using  tew  computer  invokes  a  system  error  message,  with 
cessation  of  program  operation,  so  this  condition  oast  be  recognised  and  avoided.  This 
is  done  by  checking  that  1  is  zero  and  it  is  necessary  to  use  uouble  precision  to 
ensure  that  an  exact  zero  is  obtained.  Such  cases  are  detected  by  LOCGCT,  and  Ct  is 
set  to  zero.  In  pra  tire  a  unitors  chip  is  never  created  and  none  exists  in  the  chtp 
library,  but  uni  torn  portion*  ot  images  have  been,  found  -  these  have  occurred  in  cases 
ot  scir.-Une  awthunis*  failure  when  part  ot  the  data  has  been  lost  and  the  system- 
corrected  data  has  had  a  constant  value  inserted  instead. 

the  calculation  c:  5,1  and  Ri  req-ires  a  significant  amount  of  computer  effort, 

as  *a-r.  sunnatiou  involves  >c  1  values,  each  of  which  is  the  result  of  a  multiplication  of 

two  .lumbers.  Further,  the  expression,  for  CC  has  to  be  performed  many  times  for  each 

gep  (see  he  lev'.  In  rat,  ct-.  largest  single  activity  performed  by  the  program  GCV.FIND 

is  evaluating  CC,  Ary  zwans  of  reducing  this  work  would  therefore  be  welcome.  Barnes 

ant  Sitverssr.  discuss  a  class  it  algorithms  for  fast  digital  image  registration.  It 

seems,  however,  that  their  rrethoc  is  suita&Le  for  cases  vt.*.;*  the  images  being  compared 

have  approximate’ iy  similar  brightuss,  :c  cor responding  pixels  have  similar  value.  In 

the  caw  of  Landsat  iaicery  the  pixel  values  nay  vary  considerably  from  inage  to  image 

although  the  general  pattern  returns  the  same.  and  the  fast  registration  method  does  not 

.f 

appear  tc  be  applicable.  It  is  of  interest  to  note  that  Ortt  c  :  c .  chose  to  use  the 
cross-cor.-elatior.  coetticient  ac thou  defined  by  expression  tl)  for  notching  purposes. 

The  aimer ical  value  of  a  correlation  coefficient  indicates  tew  aax> unt  of 
correlation  between  the  two  items  being  compared.  CC  values  range  from  ♦!,  indicating 
perfect  correlation,  or  identical  patterns,  to  -I,  occurring  if  one  pattern  is  the 
’negative’  of  th*  other.  CC  values  near  to  zero  ia^Iy  little  correlation.  Experience 
of  CC  values  found  with  Land  it  JtSS  imagery  is  mentioaed  in  section  -.5. 

3.1  The  correlation  coefficient  surface 

For  the  purposes  of  this  and  succeeding  sections,  it  is  useful  to  consider  the 
concept  of  a  ’surtace’  of  CC  values.  For  any  given  gep  chip,  it  is  possible  to 
calculate  CC  centred  on  all  colunxi  and  row  locations  on  the  image.  The  values  can  be 
displayed  in  three-dimensional  fora,  with  the  two  ^Horizontal  axes  being  the  coluvt  and 
row  directions  of  the  image  and  vertical  direction  being  the  value  of  CC.  This  will  be 
referred  to  here  as  a  CC  surface. 


In  practice  it  is  found,  as  might  intuitively  be  expected,  that  this  surface  has 
the  shape  of  an  irregular  hill  standing  out  from  a  surrounding  plain  having  mean  value 
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tera  but  nth  local  ooiic  leading  to  individual  racdsa  values  in  the  r 
-0.25  to  *0. 25 .  The  peak  of  tbe  hill  corresponds  with  the  position  of 
this  is  the  best  relocation  place  for  that  chip  an  that  image. 


In  practice,  the  CC  surface  may  show  a  lumber  of  other  peaks,  these  being  dam  to 
other  ground  features  having  a  similar  general  appearance.  For  example,  a  gep  centred 
on  a  small  body  of  water  such  as  a  reservoir,  will  correlate  with  other  such  bodies  of 
water  of  approximately  similar  shape  and  size.  It  is  possible  therefore,  to  obtain 
apparently  good  relocations  which  are  in  fact  incorrect. 


3.4  Belocatii 


control  point 


In  operation,  the  exact  location  of  a  chip  in  the  new  intge  is  not  of  coarse  knows 
(otherwise  there  would  he  no  need  to  find  it),  hut,  as  will  be  shown  in  section  5,  tbe 
approximate  location  trill  always  be  known,  together  with  some  measure  of  the  probable 
naxionm  distance  away  from  the  actual  location,  say  I  pixels.  These  values  - 
approximate  column  and  row  nuaber,  and  largest  expected  distance  from  the  actual 
location  -  are  passed  to  subroutine  LOOCO  as  parameters  and  used  by  it  in  its  search. 
LOCGCP  is  described  in  detail  in  section  ?.  It  operates  by  starting  at  the  supplied 
approximate  location  and  conducting  an  outward  rectilinear  spiral  search  from  there. 

The  movements  -  upward,  rightward,  downward,  leftward  and  so  on  -  are  due  to  the 
rectangular  grid  nature  of  the  data.  IDOGCF  stops  the  search  either  when  it  has  posi¬ 
tively  located  the  gep,  (as  explained  below) ,  .~t  when  it  has  completed  the  specified 
oudber  of  "turns*  of  the  spiral. 

Due  to  the  vthod  of  opcra'on  of  the  m  in  program,  it  is  acceptable  to  have  am 
occasional  mislocatiou.  Experience  with  cany  hu.  reds  of  geps  in  vctIous  image*  of 
Britain  ahows  that  the  following  algorithm  gives  an  acceptably  high  proportion  (w>re  than 
9DZ)  of  correctly  located  geps  at  the  first  attempt. 

As  the  spiral  search  proceeds  from  its  start,  a  correlation  coefficient  is 
calculated  for  each  position.  The  value  of  the  largest  CC  is  stored  (in  fact  tbe  10 
largest  values  are  stored  is  descending  order,  for  later  printout)  and  each  time  a  high**- 
CC  value  is  calculated,  the  previous  highest  value  is  replaced.  If  at  may  time  tbe 
mavimum  CC  is  greater  than  0.5  and  two  coaplete  rings  of  the  spiral  have  been  performed 
without  a  higher  CC  value  being  found,  then  the  subroutine  is  brought  to  *  halt,  at 
consideration  of  tbe  CC  surface  shows  that  it  is  likely  that  the  correct  location  has 
been  found,  and  is  new  being  passed.  In  the  absence  of  such  a  positive  relocation, 

LOCGC?  continues  isitil  all  K  rings  have  been  searched.  It  then  returns  the  highest  CC 
value  found,  provided  that  it  is  greater  than  0.3,  together  with  its  position.  If  no  CC 
higher  than  0.3  is  aet,  an  indication  is  provided  that  no  relocation  point  can  be  found. 

4  SUPERVISED  RELOCATION  OF  GRDCKD  CORTtCL  POUTS 


be  found. 


Before  tbe  operation  of  tbe  program  GCP.FIKD  is  described  in  detail  it  is  of 
interest  to  consider  the  steps  leading  to  its  present  design. 

Work  on  the  problem  of  gep  relocation  started  with  an  informal  feasibility  study. 
It  soon  became  apparent  that  a  chip  size  of  about  20  pixels  square  could  be  correlated 


It  soon  be 
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with  an  image  in  an  acceptably  abort  lo^attr  tint,  and  that  a  apxral  search  was  able  to 
relocate  geps  in  a  high  proportion  of  cases.  The  operation  was  Coras <1  into  a  counter 
prograa,  called  GO. LOCATE,  (not  described  in  this  Report)  whose  operation  was  the 
basis  for  the  present  subroutine  LOCCCP,  though  lacking  scat  of  the  feature!  now  possessed 
by  that  subroutine. 

This  relocating  prograa  vn  then  applied  to  ranges,  in  a  user— supervised  wanner. 

The  user  started  by  asking  an  estimate  of  the  location  of  the  first  gep  within  the  new 
image,  and  then  osed  the  prograa  to  perfoca  the  enact  relocation.  Another  nearby  gep  waa 
then  chosen.  Referring  to  the  original  image,  vhe  second  gep  was  at  a  known  distance,  in 
coluans  and  rows,  froa  the  first  gep.  These  colora  and  row  differences  could  then  he 
applied  to  the  location  of  the  first  gep  in  the  new  iaage.  to  give  an  approxiante  location 
for  the  second  gep  in  that  iaage.  GCF.  LOCATE  was  then  used  to  perform  the  exact  re¬ 
location.  This  operation  assumes  that  the  new  iaage  is  to  the  rat  scale,  and  is  orientated 
at  the  sane  angle,  as  the  first  iaage,  an  assuaption  which,  whilst  seldom  exactly  true . 
is  never  far  from  the  case.  The  differences  of  scale  and  orientation  clearly  introduce 
less  absolute  error  for  closely  spaced  geps,  hence  the  choice  of  near  neighbours. 

When  several  geps  have  been  located,  re  by  one,  with  the  aid  of  GCP. LOCATE,  a 
trans  format  ion  matrix  was  then  calculated,  using  least-squares  fitting.  Inspection  of 
the  residuals  -  the  distance  of  each  gep  relocation  point  from  the  ’best*  fit  given  by 
the  least-squares  calculation  -  showed  whether  any  geps  had  been  badly  aislocxted, 
(assuming  that  moat  were  correct).  Such  had  locations,  if  any,  were  then  rejected  and  a 
usv  transformation  matrix  calculated.  The  resulting  matrix  allowed  estimates  to  he  made 
of  the  expected  positions  of  other  geps,  and  these  were  then  relocated,  one  at  a  time, 
using  GCP. LOCATE.  After  several  more  geps  bad  been  located,  another  transformation  matrix 
was  calculated,  this  being  rather  wore  exact  than  the  first  one,  and  the  process  repeated 
until  all  geps  within  the  new  inage  had  been  examined. 

At  the  end  of  these  operations,  moat  of  the  geps  had  been  correctly  relocated,  but 
a  few  had  not  been.  Using  the  last-calculated  transformation  matrix,  the  expected 
positions  of  these  few  were  recalculated  one  by  one  and  new  searches  made  for  each.  This 
sometimes  resulted  in  satisfactory  relocations  being  achieved,  due  to  the  better  matrix, 
thus  adding  to  the  list  of  gep  locations. 

A  first-order  matrix  is  adequate  for  determining  the  approximate  location  of  a  gep 
and  an  existing  program  MATRIX  5  waa  used. 

When  several  images  had  been  processed  by  the  net  bod  just  described,  it  becaae 
apparent  that  the  user  i as  employing  a  ’human  algorithm*  and  chat  with  experience  the 
’decisions'  involved  were  seen  to  be  of  a  nature  amenable  to  coa^utation.  It  man  there¬ 
fore  decided  to  combine  the  sequence  of  operations,  am  far  as  possible,  into  one  all— 
embracing  pt  mm  Experience  showed  that  not  only  was  this  possible  hot  that  it  worked 
well  for  m>st  images  exrained.  With  farther  experience,  refinements  were  added  wa 
the  present  prograa  GGP.FEQ)  had  beau  evolved. 
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AUTOMATIC  KLOCATIO*  OF  CgOUXD  OOXHOL  FOOTS  IT  MPCMM  CCF.rOP 


Following  the  work  described  in  section  4,  the  prog  raw  QCF.FIMD  was  evolved.  When 
this  program  is  ran,  it  asks  the  user  a  few  questions.  Once  these  hose  been  answered, 
the  progran  then  conducts  its  search,  totonatically  relocating  the  Majority  of  gcps  and 
providing  a  list  of  their  locations.  This  progran  will  now  be  described  in  sene  detail, 
the  subroutines  being  described  in  section  7. 

5.1  Peer-supplied  infomatioo 

The  information  requested  fron  the  user  includes: 


the  filensae  of  the 
in  sane  cir 
the  identity  .t  the 
soar  indication  of 


details  shout  this 
chosen  first  gep, 

location  of  that  gep  in  the 


The  first  three  of  these  four  iteas  are  of  an  leidrnnwf i ng  nature  and  ;  Jeed  the 
second  iten  is  no  longer  requested  in  the  majority  of  cases.  Only  the  foutti  item 
requires  sane  effort  fron  the  user,  and  the  possibility  of  performing  that  automatical ly 
is  discussed  in  section  S.7.  Apart  fron  iten  4.  the  operation  of  the  progran  can  be 
considered  to  be  entirely  automatic. 


The  second  iten  listed  above  concerns  details  about  the  inage  nenrd  in  the  first 
item.  This  image  (which  nay  in  fact  be  a  subinage  of  an  entire  Landsnt  MSS  image)  ban 
associated  with  it  certain  information,  such  the  the  Inedeer  track  and  frame  amber,  the 
date  of  the  Lssdsat  pass,  and  so  forth,  llben  this  work  started,  these  details  were  not 
provided  with  the  iunge,  and  had  to  be  supplied  by  the  user  at  the  terminal  keyboard. 

It  was  therefore  decided  to  attach  a  ’tail’  to  the  end  of  each  image  file,  containing 
such  information  about  the  image.  This  tail  is  now  automatically  provided  whenever  an 
inagr  is  constructed  on  disc  storage  fron  the  cauputer— compatible  sagnetic  tape,  by  the 
program  tH.FPC.NSS*.  If  a  subinage  is  formed  fron  an  entire  image,  using  the  program 
H.SCT*  the  current  version  of  chat  progran  ensures  that  a  suitable  tail  is  attached 
to  the  subinage,  providing  in  addition  to  the  details  mentioned  above,  the  correct  top 
left  hand  corner  colon  and  row  aadert  referred  to  the  entire  inagr .  Since  inagr*  are 
nor.  usually  retained  for  long  periods  of  time  on  disc,  in  practice  meet  images  now  have 
tails,  so  the  second  item  is  aeldon  requested. 

The  choice  of  the  first  gep  to  be  located  is  at  the  user's  discretion  and  should 
preferably  be  one  which  stands  out  clearly  fron  its  surrowdiags.  It  is  also  preferable 
to  cbooae  a  gep  well  away  from  the  edges  of  the  inage  and  ss  far  as  possible  eun  i  loaded 
by  other  geps,  so  that  there  are  oeuy  close  ones,  since  this  reduces  the  searching  tine. 
In  practice  it  is  helpful  to  rhooae  s  second  'initial  gep*  at  the  sane  tine  as  the  first 
one  is  chosen,  so  that  if  the  prine  one  cannot  he  correctly  located  a  second  one  is 
immediate  ly  available. 


If  only  cloud-free  images  were  to  be  eiminrt,  one  gep  could  bn  selected  as  the 
'standard'  first  gep  for  each  track/ frane  in  the  Landsnt  system.  However  this  progran 
nay  successfully  be  used  on  ineges  which  have  considerable  cloud  cover,  provided  that  n 
proportion  of  the  gepe  are  visible,  ad  in  sack  cases  oner  discretion  in  the  choice  of  gep 
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The  user  is  asked  for  Mac  indication  M  the  position  of  Che  first  gcp.  as  the 
search  is  only  conducted  for  40  'rings'.  t£  up  to  40  pixels  in  all  directions.  The 
position  nap  be  provided  in  one  of  two  nodes,  either  as  approximate  celon  and  rear 
nosers  in  the  new  image  or  from  Measurements  taken  f  ran  a  picture  of  the  Land  sat  scene, 
the  four  dimensions  required  being: 

x  «  horizontal  distance  of  gcp  from  left  hand  side  of  picture, 

X  »  horizontal  width  of  picture, 
y  »  vertical  distance  of  gcp  from  top  of  picture, 

T  «  vertical  height  of  picture. 

Any  unit  of  distance  nap  be  used,  provided  that  the  sene  unit  is  oaed  for  nil  four 
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measurements.  In  the  case  of  a  standard  ME  Li  noscan  print,  X  and  T  are  each  a 

little  more  than  180  ms.  Provided  that  all  of  the  four  measurements  are  accnrate  to  the 
nearest  me  (:V  their  error  does  not  exceed  0.5  an)  then  the  error  in  the  implied  lo<  tins 
on  the  image  should  not  exceed  18  coltams  or  13  rows,  which  is  well  within  the  search 
spiral  of  40  rings. 

5.2  Estimation  of  initial  transformation  matrix 

When  progran  GCP. LOCATE  was  used  manual lp.  the  scale  and  orientation  of  the  aen 
image  were  obtained,  bp  inplication,  bp  assuming  them  to  be  the  same  an  for  the  original 
image.  In  the  present  program,  the  scale  and  orientation  are  calculated  from  knowledge 
of  the  Landsat  sp stem  mid  the  track  and  frame  maker.  Thus,  an  initial  tr  ana  formation 
matrix  can  be  calculated,  to  determine  the  expected  positions  of  the  second  and  a  aha  eqocat 
g  T>»*  once  the  first  gcp  has  been  located.  It  is  recognised  mat  this  matrix  is  ouly  an 
approxination,  but  it  5s  satisfactory  if  it  is  good  enough  to  enable  several  nearby  gcp* 
to  be  located. 


To  calculate  the  initial  matrix,  the  following  operations  are  done.  The  pixel  siae 
is  assumed  to  be  57  bp  79  net  res,  these  being  typical  values.  The  orientation  of  the 
imagr  with  respect  to  the  map  coordinate  system  l THETA)  is  calculated  for  that  »”gr  track 
and  frame  nunber,  bp  subroutine  U9GLE,  di scribed  in  section  7.  The  rotation  of  both 
x  and  p  direction  are  taken  to  be  of  value  TKTA,  which  inplies  that  the  image 
is  not  sheared.  The  angle  of  shear  depends  on  the  paw  angle  of  the  satellite,  i*  map 
deviation  of  the  sea  direction  from  perpendicularity  with  the  direction  of  satellite 
travel,  conbined  with  earth  rotation  effects,  and  it  has  proved  satisf actorp  to  assume 
zero  shear.  Using  these  aseunptions,  four  of  the  six  nstrix  values  can  he  calculated  and 
placed  into  arrays  81  and  82.  These  four  values  determine  the  scale  and  orientation  of 
the  imegc  in  relation  to  the  sap  coordinate  system.  The  final  two  aatiix  values  determine 
the  relative  displacement  and  are  provided  by  the  operation  described  in  section  5.3. 


5.1  lelscatiou  of  first  ground  control  point 


As  the  user-supplied  information  is  provided,  severs!  starting  operations  are 
performed  bp  the  program,  as  follows. 


Once  the  filename  for  the  new  image  has  been  supplied,  that  file  is  opened. 

open  correctly,  eg  a  non-existent  filename  has  bean  supplied,  invokes  a  rupect 


Failure  to 
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request  for  a  aac.  Once  opened,  the  file  is  examined  to  see  whether  it  has  a  tail.  If 
it  has,  relevant  data  is  extracted  (ran  it,  it  not  the  information  is  requested  from  the 
user  —  Laid  sat  track  and  frame  number  and  pass  dace. 

After  the  required  first  chip  identity  has  been  Supplied,  in  the  form  of  mnp  mi 
serial  nuabers,  subroutine  CBFSASf  is  called,  to  construct  the  filename  for  the  gcp  chip. 
An  invalid  chip  identity  is  recognised  by  CBFSdH.  and  a  repeat  request  is  made  until  a 
satisfactory  identity  is  provided.  The  master  map  location  file  is  them  searched  for 
this  chip  identity  and,  if  found,  the  rr lever'  .iap  easting  and  northing  are  stored  and 
the  chip  file  opened.  Failure  to  locate  the  .nip  in  the  master  nap  file  camel  •  request 
for  another  chip  identity. 


An  attenpt  is  then  nade  to  relocate  the  first  gcp  chip  specified .  If  the  approxi¬ 
mate  location  has  been  provided  in  the  form  of  four  picture  dimensions,  (section  S.l)  then 
the  program  calculates  the  approximate  coluun  and  row  numbers.  This  calculation  ss eases 
that  the  pictorial  representation  of  the  image  is  in  the  form  no  reel  ly  available  in 
Space  Department.  U£,  :V  a  'quick  loot’  or  better  a  Linoscan  print,  ir.  which  esses  the 
whole  OCT  image  is  used,  resulting  in  a  slightly  skewed  parallelogram.  The  asemt  of 
skew  is  approximately  one  pixel  in  twenty  for  most  images  of  Britain,  with  system 
corrected  data.  If  the  image  is  not  system  corrected,  or  if  the  image  is  at  a  different 
earth  latitude,  then  the  amount  of  skew  will  probably  differ.  In  such  esses  the 
calculation  of  the  column  number  would  need  modification. 


As  an  alternative  to  providing  the  four  picture  dimensions,  the  user 
approximate  location  as  Culiss  and  row  raises,  and  mo  conversion  is  needed. 


supply  an 


Subroutine  L0OCCP  is  then  called,  to  request  a  search  from  the  location  provided, 
for  up  to  40  rings,  an  acticv  indicator  being  set  to  ?.0  for  a  nearest-integer  search. 
If  this  first  gcp  is  uot  found  by  LOGGGP,  indicated  by  the  parameter  KSCII  ''tuned 
less  than  zero.  message  is  displayed  at  the  user  terminal,  mod  the  program  returns  to 
an  earlier  stage,  to  request  another  gcp  and  location.  The  osar  may,  if  be  so  wishes, 
supply  the  same  gcp  but  with  a  different  starting  location,  if  this  seams  desirable. 


Once  the  first  gcp  has  been  relocated,  its  coins  md  row  numbers  can  be  need, 
together  with  its  nap  easting  and  northing,  to  calculate  the  final  two  values  of  the 
initial  first-order  transformation  matrix,  stored  in  arrays  B1  and  B2. 

5,4  Be  location  of  other  ground  control  points 

After  the  first  gcp  has  been  successfully  located  and  the  initial  transformation 
matrix  prepared,  the  relocation  of  other  geps  can  camuence.  Since  the  transformation 
matrix  is  known  to  be  only  an  approximation,  the  accuracy  of  prediction  of  the  position 
of  other  geps  will  become  progressively  worse  as  the  distance  of  the  gcp  from  the  first 
one  increases.  TSe  next  operation,  therefore,  is  a  scan  through  the  master  map  location 
file  to  list  all  cp *  which  seen  likely,  on  the  besis  of  the  initial  matrix,  to  he  inside 
or  near  tc  the  inage  lor  subinage)  in  use.  This  ’batch’  of  geps  is  tabulated  in  a  set  of 
arrays,  which  bold  the  nap  and  serial  m^ers,  eastings  and  northings  and  distance  from 
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the  first  gep ,  other  table  values  being  set  to  aero  for  later  use.  llbon  this  operation 
is  cnpletc,  the  user  is  informed  of  the  maker  of  geps  in  the  hatch,  typically  SO  to  >00. 

Each  gep  in  the  hatch  is  then  ecamimeo,  is  ascending  order  of  distance,  is  a 
similar  manner  to  the  first  one.  It*  filename  is  constructed  by  subrout ine  QfMli,  the 
approximate  location  is  calculated,  non  using  the  easting  and  northing  together  with  the 
initial  transformation  matrix,  the  chip  file  is  opened  and  L0CCCP  is  called  to  relocate 
it.  Since  the  estimated  position  of  each  gep  should  now  be  knows  much  more  accurately 
chan  the  first  gep,  a  much  smallar  spiral  search  is  sufficient,  the  maker  of  rings 
specified  being  dependent  epon  tie  distance  of  the  gep  from  the  initial  one.  In  practice 
this  is  found  to  he  Suitable,  provided  chat  the  first  gep  has  bean  correctly  located: 
if  it  has  not,  then  there  is  little  likelihood  of  any  others  being  correctly  f.<  «d.  When, 
as  usually  occurs,  the  gep  has  been  successfully  relocated,  lOOIXf  provides  details  of 
the  position  and  correlation  coefficient ,  CC,  and  the  main  program  informs  the  user  that 
the  chip  has  been  found.  It  may  be  noted  in  pasting  that  a  considerable  amount  of 
information  on  the  progress  of  the  program  is  provided  at  the  osar's  terminal :  this  may 
be  ignored,  or  a  general  watch  kept  on  the  progress,  or  the  details  can  be  saved  on  a 
logging  file  for  subsequent  detailed  inspection.  A  later  version  of  COP. F ISO  omits  mack 
of  this  detail,  and  condenses  the  remainder. 

After  10  geps  have  'keen  relocated,  with  CC  of  0.6  or  more,  the  subroutine  HkTSUB 
is  called.  This  subroutine  calculates  a  new  first-order  transformation  matrix,  based  am 
the  actual  positions  of  the  10  geps  and  this  should  represent  an  iaprovement  upon  the 
initial  transof matron  matrix,  which  was  not  based  on  actual  gepe.  Following  this 
calculation,  subsequent  UOCCCr  calls  only  <u  for  12  riogs  to  he  searched,  as  the  better 
transformation  matrix  now  available  increases  the  accuracy  of  the  estimated  positions. 

All  the  remaining  geps  in  the  hatch  are  examined,  with  a  12-ring  search,  calculating 
a  fresh  trans  format  ion  matrix  after  every  tern  more  have  hewn  relocated.  When  all  in  the 
batch  have  been  examined,  another  transformation  matrix  is  calculated.  The  Subroutine 
HUSCt  is  described  in  section  7.3,  but  it  may  he  mentioned  here  that  a  limit  to  the 
accepted  residual  size  may  be  sf~  .  .  ^  a  paraarter  to  this  subroutine,  and  thus  the 
occasional  mis  located  gep  is  not  al  lowed  to  falsify  the  transformation  matrix. 

5.5  Ke checking  of  doubtful  locations 

To  complete  the  initial  relocation  process,  poorly  located  geps  are  rechecked.  The 
subroutine  NATSCB,  apart  from  calculating  a  transformat  ion  matrix,  assigns  to  each  gep  an 
'error*  which  is  the  residual  resulting  from  that  calculation.  Carnally  soon  of  the  gepa 
have  residuals  which  are  rather  larger  than  might  be  expected,  this  limit  being 
arbitrarily  specified  as  3  pixels.  Such  geps  are  therefore  rechecked,  by  a  repeat  of 
the  relocation  process.  It  is  sometimes  fond  that  a  gep  can  ha  relocated  at  a  position 
more  consistent  with  the  others,  usually  because  the  previous  search  for  that  gep  started 
too  far  away  to  find  it.  or  perhaps  because  a  lower  peak  in  the  CC  surface  had  bean  met 
previously.  The  smker  of  such  batter  relocations  is  usual ly  small,  perhaps  two  per  ^ 

image,  dependent  on  the  quality  of  that  image.  Am  image  with  hate  is  amongst  the  meat 

esoam,  leading  to  many  mis  locations  which  have  apparently  acceptably  high  CC  valtna. 


trouble 


viu  the  existing  relocation  position  is  wtvij,  and  liXOCT  is  called,  this  tine  with  a 
limit  of  ten  rings,  since  the  position  is  now  known  quite  accurately.  The  action 
indicator  for  LOCCCP  is  set  to  1.0,  to  indicate  that  the  search  is  to  0.1  pixel.  In  fact 
LOCCCP  performs  two  spiral  searches  in  this  case-,  first  a  (maxim*  of)  six  ring  search 
for  the  integer  location,  followed  by  a  fractional  search  of  up  to  ten  rings.  Occasionally 
the  location  is  modified  during  the  integer  search,  and  this  is  usually  for  th  better  as 
judged  by  the  resuiting  residual. 

Finally  the  program  outputs  its  results  in  a  ’permanent*  form,  M  other  than  the 
transient  output  to  the  terminal.  The  final  relocation  positions,  to  the  nearest 
0.1  pixel,  are  output  in  a  *S1D  list*,  one  position  for  each  map  identity.  Additionally, 
a  final  transformation  matrix  is  calculated  and  output  as  a  file,  hut  this  is  not  often 
used  as  such,  and  has  been  omitted  Iron  subsequent  versions  of  the  program.  In  practice, 
users  tend  to  prefer  to  accept  the  StD  list  and  calculate  a  transformation  matrix  of 
either  first  or  higher  order  under  their  om  control  using  a  program  OCPFXT  Coot  described 
here). 

He  program  ends  by  closing  all  files.  The  user  is  then  asked  whether  he  wishes  to 
use  the  program  again. 

5.7  Possible  future  developments 

This  section  mentions  possible  future  developments  of  the  program  CCF.FHD.  Since 
smr  of  these  dewlopnrnts  do  not  appear  to  be  promising  from  an  overall  cost- 
effectiveness  viewpoint,  mad  use  in  other  regions  is  not  at  present  required,  no  work 
is  at  present  taking  place.  Some  of  these  ideas  night  be  useful  however,  if  die  scope 
of  the  work  were  to  be  extended,  to  cover  substantial  latd  areas  oitiidc  Britain, 
particularly  places  having  little  cloud  cover. 

5,7.1  automatic  relocation  of  the  initial  groimd  control  point 

Each  Landsat  image  has  its  centre  ai  a  defined  nominal  location  on  the  earth’s 
surface,  as  described  in  section  7.2.I(ii).  For  reasons  described  in  Bef  2  in  practice 
the  centre  of  an  image  can  be  tax  distance  from  the  nominal  position,  by  op  to  about 
50km.  It  would  be  possible  to  relocate  the  initial  gep  by  allowing  it  to  perform  a 
spiral  search  sufficiently  large  to  cover  30  km  in  all  directions  from  the  expected 
position.  This  would  require  about  300  rings  in  the  along- track  direction  and  ABO  in 
the  cro  —track  direction.  A  search  as  vide  as  480  rings  i*  open  to  two  objections,  the 
mrt  obvious  one  being  the  tine  required:  this  might  he  as  long  as  8  hours  for  the 
entire  search,  although  often  the  gep  would  be  found  before  all  of  the  rings  had  been 
examined  The  other  objection  is  based  on  relocation  experience,  which  indicates  that 
a  im^i r i  of  possible  relocations  are  likely  to  be  found  in  such  an  extensive  search,  and 
it  would  often  be  difficult  to  determine  by  machine  which  is  the  ’correct*  one  -  peak  CC 
value  alone  not  necessarily  being  sufficient. 
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A  possible  method  would  be  to  conduct  p  re  l  ini  nary  wircbti  at  a  different  'scale*. 
For  ezaaple,  an  area  190  by  190  pixels  in  site  could  be  reduced  to  a  19  by  19  chip  by, 
for  enable.  averaging  subsets  of  10  by  10  piaels,  or  by  selecting  every  tenth  pixel  of 
every  tenth  row.  The  resulting  'superchip*  could  then  be  relocated  by  a  spiral  search, 
cowparing  it  with  aiailarly  treated  portiooa  of  the  inage,  eaN  turn  of  the  spiral  being 
ten  pixels  in  width.  It  would  be  necessary  to  select  the  superchip  so  that  it  contained 
a  distinctive  pattern  appropriate  to  its  aise.  Once  this  prelimir.-ry  relocation  had 
been  wnde.  a  normal  relocation  process  could  then  be  conducted  to  obtain  the  required 
accuracy.  If  necessary,  a  three  or  four  stage  operation  cow'd  be  conducted,  using 
several  types  of  s’jperchip,  with  reduction  factors  of  say  7  >  md  5,  though  it  seena 
likely  that  awe  size  only  would  suffice. 

The  system  would  therefore  involve  the  use  of  one  superchip,  or  possibly  a  aeries 
of  two  or  three,  for  each  land  sat  scene.  This  system  would  seem  well  suited  to  regions 
of  the  world  which  Lave  little  cloud  or  base.  i!a  Britain  there  is  considerable  cloud 
cover,  and  apart  from  co^letely  musablt  scenes,  ther»;  are  many  where  a  useful  aamoat 
of  land  is  visible,  but  often  not  (say)  the  centre.  It  would  therefore  be  necessary  to 
have  a  umber  of  superchips  for  each  scene,  and  the  user  could  select  the  anst  premising 
one.  Even  so.  the  larger  scyerchips  may  no:  be  relocatable,  when  the  conventional  siaed 
ones  cm  still  be  relocated  in  the  spaces  between  clouds. 

So  experience  has  yet  been  gained  with  the  use  of  superchips.  It  does  seem  possible 
however,  that  there  may  be  difficulties  in  selecting  suitably  distinctive  areas  in  some 
inland  parts  of  Britain  and  in  other  countries,  though  the  extensively  convoluted  coast** 
line  of  Britain  anemia  provide  many  suitable  shapes. 

In  susuary  therefore,  this  method  should  be  suitable  for  some,  but  not  necessarily 
all,  parts  of  Britain  and  other  countries. 

5.7.2  Pee  of  raw  data  for  first  amd  relocation 

All  imagery  used  for  this  work  baa  been  in  the  system-corrected  font*,  which  is 
known  to  include  several  types  of  alteration.  For  images  of  Britain  the  alterations 
include  both  a  displacement  of  lines  by  one  pixel  at  intervals  of  about  20  scan  lines 
(the  exact  atmdter  varies),  and  the  insertion  of  extra  pixels  at  intervals  which  nay  he 
as  close  as  one  per  29  pixels  in  the  centre  of  each  scan  line.  Ttaos  moat  chips,  being 
19  by  19  pixels  in  size,  will  include  at  least  one  of  these  types  of  alteration.  These 
alterations  are  valuable  for  providing  am  overall  rectification  of  tbs  image,  hat  are 
made  at  the  expense  of  local  discontinuities.  It  seem  probable  lhat  both  the  initial 
location  of  geps  oo  naps  and  jlm  the  mb  sequent  automatic  relocation  could  be  done 
with  a  higher  level  of  correlation  if  such  discontinuities  did  not  exist. 

There  is  now  s  considerable  investment  of  effort  in  the  existing  set  of  chips 
extracted  from  system-corrected  images,  together  with  their  colom  sad  row  details,  so 
it  is  mot  economic  to  change  now  to  the  use  of  rsw  data.  However,  if  a  asm  system  mere 
being  set  #,  either  for  other  regions  of  the  world  or  for  a  different  type  of  sensor, 
other  chan  land  sat  HSS,  serious  consideration  should  be  given  to  the  wee  of  data  contain- 
ing  as  little  rectification  as  possible.  Since  the  nature  of  a  gep  system  is  to  allom 
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bn)mt  at  tte 


the  application  of  soar  for*  of  traafomtion,  using  ran  data  maid  wurely  rrnlt  in  a 
si ightly  wrt  coapliated  transformation  over  wbicl  oat  has  complete  control ,  ratter 
transform^  natcrial  which  has  already  suffered  ooc  prccses  of  rectification. 

■scent  informtiou  indicates  that  so*e  system  correction  is  non  teiag  done  by  a 
reselling  tecteique,  ratter  than  by  addins  pixel* .  Tte  shorn  mentioned  objection  mould 
be  leas  valid  in  this  csss. 

5.7.3  Wole  Hint  of  satallits  teteeionr 

A  conplete  rectification  systen  would  take  account  of  tte  wwral  gaoma tribal 
considerations  normally  used  for  rectification,  and  tall  also  t^e  account  of  tte 
■ncatati  of  tte  satellite,  in  particular  its  pitch,  roll  and  yaw.  Shliea  tea  stem 
that  it  is  possible  to  deduce  evidence  ah  ant  tte  satellite  pitch,  roll  and  yam  from  tte 
details  of  an  inage  if  sufficient  geps  are  selected.  Conversely,  a  nore  accurate  fit 
of  gepo  to  should  be  obtainable  if  information  were  available  stent  tte  satellite 

■irranrnrr  In  fact  this  inforaatioa  is  tenuu,  if  not  to  s  high  level  of  accuracy.  but 
the  infomation  is  not  transferred  to  CCT  (at  least  those  available  f 
tine  of  writing). 

To  obtain  the  highest  accuracy  possible,  it  uould  therefore  sewn  desirable  to 
raw  (u  not  systen  corrected)  data,  and  to  bane  available  and  nan  all  team  utfoii 
on  tte  satellite  attitude  and  ita  changes. 

5.7.4  Autowatic  decision  of  relocation  quality 

Whilst  ibserving  operation  c  tte  pragma  OCF.IIID,  tint  user  finds  hi  use  If 
judgeaEotx  about  tte  quality  of  each  relecntic-  position.  These  judgments  ere 
such  factors  as  the  value  of  tee  correlation  coefficient  CC  and  tew  clone  relocation 
is  to  tte  expected  place.  For  tinplt,  a  CC  value  of  0.9  within  one  or  cm  purls  of 
its  expected  location  appears  'good*,  whereas  a  siadlarly  high  CC  value  of  nsy  10  pixels 
nwar  Iron  its  expected  place  excites  soar  doubt.  Bonnily  it  is  tetter  to  suspend 
judgemot  until  tte  whole  iwagt  has  teen  processed,  st  which  stage  the  poorer  locations 
are  re checked. 

Tte  author  considers  that  it  is  tetter  to  restrict  tte  operation  of  tte  progran 
GCP.FI5D  to  actually  finding  tte  locations,  sad  leave  the  jodgewnat  of  their  quality  to 
a  separate  progran,  (such  as  GCTFIT)  applied  to  then.  Tte  only  occasion  during  progran 
operation  when  a  judgment  is  needed  is  daring  tte  early  stages  thn  n  ratter  higher 
value  of  CC  (0.6)  is  insisted  upon,  to  reduce  tte  possibility  of  iacorract  locations 
upsetting  tte  first  aatrix  recalculation. 

5.7.5  legions  of  tte  world  otter  than  Brit  sin 

CCF.FIXD  was  designed  originally  for  use  with  Lsndsst  inages  of  Britain.  For 
exaaple,  tte  naps  used  are  in  the  British  Batioaal  Grid,  which  is  not  applicable  to 
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otter  parts  of  tte  world. 

Tte  progran  could  be 
widespread  nap  system,  such 


rsl  in  its  application  hy  using  sob 

lercator  (CM),  this  would  chiefly 
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involve  alterations  to  tbs  subroutine  Itfl5  and.  depending  iqmo  the  nap  systen  need, 
night  require  additional  data.  Iron  the  user .  such  as.  for  regions  ness  the  edge  of  a  DIM 
zone,  tbicb  zone  *s  to  he  used.  The  nap  reference  naster  file  of  gepo  would  of  course 
have  tv  he  presided  in  the  appropriate  coordinate  systen. 

A  few  other  changes  would  be  necessitated  and  these  have  been  motioned  during  the 
course  of  the  prcgvws  description. 

6  ormrm;  «  npqug  ocr.rna 

This  section  provides  detailed  instructions  on  how  to  operate  the  progran  GCP.FDO. 
It  as  tones  sons  faniliarity  with  the  boat  ca^uter  systen  CPrin  750)  sod  cane  experience 
of  using  canputer  progran  in  general. 

6.1  Initial  coeditions 

The  progran  make*  use  of  (a)  a  Landaat  H5S  inage  in  the  forn  of  n  Space  Department 
inage-fomat  file,  Cb)  a  neater  nap  location  file  and  (c)  a  set  of  gep  chips,  so  all  of 
these  neat  be  available  before  operation  can  coassmce.  I  tens  (h)  nd  (c)  are  normally 
present,  in  the  user  areas  required  by  CO.FBD,  but  it  is  the  user's  responsibility  to 
provide  item  (*) -  Normally  a  lietsat  band  7  insgrr  is  used,  because  ell  of  the  chips  have 
been  prepared  f ron  hand  7  inages.  Inages  in  other  bands  have  been  tasted,  hot.  unsurpris¬ 
ingly,  band  6  has  shown  less  successful  results  than  band  7,  and  bards  5  and  4  have 
shown  pc  or  results.  Normally.  an  image  is  obtained  in  the  forn  of  a  m^iiiisr  rnspetihli 
nagnetic  tape  (OClj  containing  all  four  hands,  so  selection  of  bud  7  presents  no  problem. 
Innge-fo*  snt  k’i.lej  any  oe  obtained  t'rui  a  CC1  by  use  r.l  the  profxan  IM.nK.lBS4,  in 
which  case  a  tail  is  provided  for  each  file,  for  purposes  of  coeparison,  it  nsy  he 
menriooed  that  the  typical  operating  tine  to  transfer  an  ienge  1/  nee*-*  of  ra.mC.N5S 
is  ot*«t  one  aour,  since  it  involves  two  passes  of  a  2400  ft  taps,  .t.jis  figure  iimnil 
use  at'  the  systen:  opei  j.tinu  siy  he  slowed  in  *  multi-user  environment. 

6.2  Ope.vicQ 

Progrns  OL'.'tSD  has  he’s  observed  to  relocate  all  geps  in  a  new  rugr  in  a  typical 
sol  -user  elapsed  tine  of  shoot  one  quarter  minute  per  gep,  w  25  somites  for  about 
100  g.-ps.  The  progrm  has  been  designed  so  that  all  user  interaction  takes  piece  st 
the  beginning,  after  which  the  process  nay  be  left  unattended.  It  is  recommended  that 
information  presented  at  the  user  terminal,  nocb  of  which  passes  rapidly,  should  he 
1  «*ed.  so  that  it  cm  he  inspected  at  leisure,  if  need  arises,  ft  nsy  he  convenient 
to  ns  the  program  as  a  'phaeton* ,  to  frse  the  tcrninal  for  other  use,  but  if  this  is 
done  it  is  desirable  to  ensure  first  that  the  initial  locations,  up  to  the  first  matrix 
recalculation,  are  satisfactorily  achieved. 

Fig  I  lists  the  first  45  lines  of  operstion  seen  on  the  terminal  daring  n  typical 
run.  This  figure  will  now  be  considered  line  by  line. 

OOMD  L0CFIL  This  user  raciest  npm«  a  log  file  named  LOCFIL  which  logs  til  of  the 
subsequent  terminal  output,  llhen  the  run  has  been  co^lcted  the  oner  should  key  in 
COMD  H9D  to  clone  the  current  log  file. 
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GCT.FUD  ibis  rtqo«(ti  the  rrogrM  to  be  no.  TV*  frogria  aaaosncea  its  nans  and  wrtki 
date.  Tbs  current  date  aud  list  an  then  output ,  together  with  two  tiwss  vbicb  rsfsr  to 
tbs  processor  and  input /output  tins  wasps. 

Dill  Wl  Tbs  initialising  of  tbs  DUB  system  by  was  of  scA  routine  IBIT1  •  section  7.$) 
caiMt  this  nessage  to  bs  output. 

nDWIDC  DETAILS  OF  VX.  SUP: 

TTPC  l*TL*T  FILE  XAME:  A  'snap*  Up  lies  a  unipus  l  anils  at  t  rack/ f  ms /data.  Tbs  reply 
should  be  tbs  news  of  tbs  inegr-fomnt  fils  containing  tbs  band  7  inapt  which  is  to 
base  relocations  performed.  If  an  invalid  or  non-existent  filename  is  provided,  a 
further  reguest  will  bs  made.  Ubsn  the  fils  has  bssn  opened,  its  "teat*  is  output, 
together  with  its  siss  in  columns  and  rows.  In  the  present  cans  the  inapt  fils  has  a 
tail,  otherwise  the  user  would  bs  asknd  to  provide  tbs  track,  frsan  and  data. 

IXITLAL  GBOCXD  COTTXOL  POUT :  XXX  X  Tbs  user  responds  with  bis  selected  initial  pep 
identity  in  tbs  form  of  two  integer*  r  the  usher  and  serial  sunhsr .  If  tbs  identity 
is  illegal  or  does  not  exist,  tbs  user  will  he  informed,  otherwise  tbs  chip  fils  will  be 
opened,  its  'text*  written  out  nd  its  size  provided  -  this  should  always  be  19  by  19, 
as  hare,  with  the  present  system. 

mkjvidc  position  of  ccr  aid  size  of  fictcxe 

DATA  (Ml.  BEAL)  n  FOB!  iXyT.  Os  IF  COL  AB  KM  Tbs  user  any  provide  picture 
dimensions  as  described  in  section  5.1  or,  as  bar*,  provide  an  ■pprnsinata  colons  and 
row  number.  In  the  latter  case: 

AFT’BDXDIATE  C.C.P.  UOCCATIOS  BETTBXF.D  TO  EXT1BE  DVCE 


OOL:  User  provides  appro vi ware  cole 


SOU:  User  provides 


eroxiaetc 


Another  version  of  the  program  rv guests  the 
list.  For  this  version,  a  standard  temporary  fils 
ntchimc. 


of  the  file  for  the  ontput  *SID" 
is  entwastically  provided  by  the 


The  program  nay  then  be  left  to  conti one,  proceeding  in  the  meaner  described 
elsewhere,  until  couplstioo.  Fig  I  shows  the  spiral  search  for  the  first  gep  sad  its 
successful  completion.  Finally: 

AXOTHEB  MMX?  to  which  the  user  replies  TES  or  XO  as  he  reguires. 

A. 3  Ei Tor  conditions 


As  already  mentioned,  operation  nay  not  he  trouble-free.  During  the  above 
description  it  was  noted  that  invalid  replies  to  questions  often  gave  rise  to  error  or 
waning  messages-  Other  possible  causes  of  error  messages  are  now  described. 

A. 3, 1  Messages  from  the  nsia  program 


The  following  aessages  cm  be  provided  by  the 
effects  are  explained- 
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FAILED  TO  OF  I*  NAPGCP.  MASTER  The  niltr  up  location  file  «u  not  ivcctMfulljr  op— 4. 
la  practice  this  should  newer  happen.  Causes  could  he  (a)  a  ijitta  fault  of  soon  kind  or 
(b)  the  Master  file  is  not  present  -  perhaps  it  has  been  deleted.  Since  the  program 
c snoot  proceed  without  the  information  contained  in  this  file,  it  stops.  The  user  is 
advised  to  consult  the  system  manage. . 

THIS  SHAT  DOES  NOT  HAVE  A  "MAIV  TAIL  This  message  will  occur  if  so  image  is  used  which 
has  been  in  existence  since  before  tails  were  invented.  In  this  case  the  user  must 
answer  some  questions  -  track,  frame  and  date  -  the  answers  to  which  are  normally 
provided  by  information  in  the  tail. 

THIS  CO*  NOT  IX  MATCCP  FILE  Occurs  if  the  user  provides  a  gep  identity  (map  and  serial 
numbers)  which  does  no.  exist  .  The  user  mist  provide  another  pair  of  values 
corresponding  to  an  existing  gcp. 

GCF  nollCI  The  prog:.'£m  makes  some  checks  on  the  size  of  each  chip,  when  it  is  opened. 
All  chips  are  19  by  l>  p>.wls,  and  should  never  fail  this  test.  The  message  has  never 
occurred  in  practice,  li  it  did  occur,  the  chip  in  question  would  he  ignored  and  a  new 
'initial  gcp*  would  be  requested. 

THIS  GCP  CAKSOT  HE  LOCATED.  TIT  AMOTSEt  If  the  initial  gcp  ia  not  relocated  by  LOCGCr 
within  its  AO-ring  search,  this  nr s sage  is  output.  The  user  should  select  another 
initial  gcp  or  provide  another  estimated  location  for  the  present  oat,  repeating  the 
process  if  necessary  until  relocation  ia  successfully  achieved. 

6.3.2  Messages  from  subroutine  LOCCCF 

MS  IS  TOO  MSMT  EISGS  TOS  LOCCCF  This  would  occur  if  lOCCC?  were  called  with  a  nahr  PP 
of  rings  in  excess  of  41.  Should  any  of  the  declared  parameter  values.  HSQ  etc,  he 
changed,  the  value  41  might  be  affected.  This  message  is  really  an  aid  for  the  progrmi 
writer,  as  it  would  expose  certain  program  writing  errors.  Users  should  never  ace 
this  message. 

KHWUUOI  (S  NOT  POSITIVE  This  message,  too,  ia  a  survival  from  the  program  writing 
stage,  and  the  message  should  now  never  be  seen  by  users.  During  the  correlation 
coefficient  calculation  the  square  root  of  14  -  12  i  13  (section  3.2)  ia  calculated. 

•A  must  not  he  negative  as  otherwise  the  system  routine  SQKT  will  fail  and  BA  must 
not  he  zero,  as  otherwise  the  resulting  CC  would  he  infinite  (the  result  of  a 
’divide  by  zero*.  This  message  informs  the  user  if  either  of  these  conditions  should 
occur. 

6.3.3  Message  from  subroutine  HATS  UP 

MATS CH:  OBLT  K  SUITABLE  GCPS  FOOCD  MATSCB  cannot  perform  a  meaningful  calculation  if 
less  than  four  geps  have  been  relocated,  so  the  program  is  brought  to  a  stop.  This 
occurrence  is  related  to  the  nature  of  the  image  or  the  position  of  the  initial  gcp 
within  it.  For  example,  if  for  sosm  reason  the  initial  gcp  has  been  mislocated,  the 
subsequent  gcp*  will  probably  be  either  mislocated  also,  or  not  relocated  at  all.  Thus 
there  will  he  insufficient  apparently  good  geps  for  HATSUB.  Again,  if  there  are  few 
clear  geps  on  the  image,  due  to  the  presence  of  clouds,  only  a  email  number  may  be 
found  in  the  batch  of  those  within  50  km  of  che  first.  In  anat  cases  the  user  can 
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deduce  the  nature  of  the  protln  (ton  inspection  of  the  log  fill  of  the  program  run  and 
e nomination  of  a  picture  of  the  image. 

6.3.4  Message  from  subroutine  flMJf 

UKU:  PATH  M  OH  KM  (I  OCT  OF  RAKE  Subroutine  flMf  checks  the  validity  of  the 
ludut  path  (track)  and  row  (frame)  within  the  limits  of  the  whole  system.  1-251  and 
I-II9  respectively.  Since  HIM F  is  designed  for  use  in  Britain  only,  it  would  perhaps 
be  wore  apposite  if  these  bounds  were  restricted  to  allow  only  British  scenes  to  be  dealt 
with.  This  message  should  never  occur  whew  using  n  image  having  a  tail,  as  the  track 
and  frame  nuabers  should  always  he  vaHd.  In  the  case  of  a  tail-lets  image  the  user  is 
asked  to  supply  these  details  and  could  provide  invalid  valuns.  in  which  case  the  message 
would  be  output  and  the  program  would  atop. 

6.3.5  Other  messages 

Host  other  messages  provide  information  on  the  operation  of  the  program  rather 
than  about  errors.  A  complete  run  of  GCP.PHD  produces  a  log  file  which  may  be  many 
hundreds  of  lines  long  and  provides  considerable  detail  concerning  the  progress  of  the 
search.  It  is  hoped  that  moat  of  these  messages  are  salt  explanatory.  As  mentioned 
above,  later  versions  of  the  program  provide  less  output  information. 

7  SUMOCTUES  USED  BY  PBOGBAM  CCP.FDTO 

7.1  Subroutine  LOCOCP 

The  general  method  by  which  subroutine  LOCOCP  determines  the  location  of  a  grp 
chip  in  a  new  image  has  bees:  described.  This  subsection  describes  the  subroutine  itself 
in  soee  detail. 

Subroutine  LOCOCP  has  four  perameters;  the  column  and  row  uunbera  at  whicb  the 
searct.  is  to  start,  the  nunber  of  rings.  R,  to  search  and  an  action  indicator.  After 
operation,  it  returns  the  relocation  tolian  and  row  nuabers,  if  found,  and  the  result 
of  the  search,  which  any  be  in  the  form  of  a  correlation  coefficient  value,  but  if  less 
than  zero  in  value,  is  interpreted  as  indicating  the  miner  in  which  the  setrch  failed. 
Thus,  if  the  best  C€  value  found  is  less  than  0.3.  this  is  not  considered  to  be  a 
relocation,  and  the  indicator  is  returned  as  -I,  to  indicate  this  fact.  A  value  of  -2 
is  returned  if  the  subroutine  is  asked  to  search  starting  from  a  location  which  is 
outside  tbs  image.,  or  insufficiently  within  it:  any  positim  less  than  9  pixels  in  fxom 
any  edge  of  the  inage  will  not  allow  c  correlation  between  chip  atafi  image  to  be 
calculated.  It  is  the  responsibility  of  the  calling  program  to  examine  the  result 
returned  by  LOCOCP,  and  to  taka  appropriate  action. 

The  next  operation  of  LOCOCP  is  to  read  the  chip  pixel  values  into  on  array 
1PA1CB,  of  sire  19  by  19.  As  each  line  of  the  image  is  rend  in,  the  opportunity  is 
tnien  of  performing  the  calculation  of  RZ  (referred  to  in  section  3.2)  for  use 
subsequently  in  the  CC  calculation. 

7-1.1  The  search  workspace 

A  Lands  a  t  }CS  image  is  normally  over  3000  by  2030  pixels  in  six*,  which  is  far 
too  large  to  hold  in  the  fast  memory  of  the  host  counter  (a  Prime  750  wscfcine),  so  the 
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data  is  bald  available  iu  alow  memory.  It  ia  therefore  advantageous  to  transfer  a 
portion  of  the  iatp  into  fast  memory  a  d  then  conduct  the  spiral  search,  rather  than  to 
perfona  the  search  by  repeated  access  oL  data  fro*  slow  enmory.  An  early  operation  of 
UOCOCT,  therefore,  is  to  transfer  a  suitable  nwat  of  the  inage  into  the  work  space, 

IHNK,  in  fast  nenory ■  The  size  of  WOK  has  been  chosen  as  101  pixels  square,  ie  a 
little  over  1 0000  words,  which  is  s  reasonable  size  for  fast  nenory .  This  square  array 
is  centred  on  the  starting  location  and  thus  allous  41  rings  of  the  spiral  to  be 
searched.  If  an  1  value  greater  than  this  ia  requested,  an  error  Message  is  output. 

If  an  1  value  less  than  41  is  specified,  only  the  required  amount  of  iusge  is  brought 
into  the  works  pare,  minimising  the  data  transfer  time.  Thus,  if  the  location  ia  known 
fairly  accurately,  the  whole  process  is  speeded  up,  both  because  of  the  reduced  amount 
of  data  to  transfer  (s  slow  input  process)  and  the  reduced  computation  tine  involved. 

If  the  specified  search  location  is  close  to  the  edge  of  the  inagr  (though  not 
of  course  f<ir*r  than  9  pixels)  there  any  not  be  a  complete  square  of  imsgi  around  that 
location.  This  situation  can  be  handled:  L0CGOP  reads  in  as  much  image  as  is  available 
into  the  works, .ice  and  then  conducts  tlae  search.  If  at  any  time  during  the  search  the 
chip  extends  beyond  the  available  image,  then  no  CC  calculation  is  made.  This  arrangement 
ensures  lust  the  maximum  .mount  of  searching  which  is  mnthemmticelly  meaningful  is  in 
fact  done. 

7. 1 .2  The  spiral  search 

Kith  both  chip  and  search  are*,  of  image  available  is  fast  memory,  the  outwardly 
spiralling  search  can  cuemnce.  The  search  starts  at  the  location  specified  by  the 
input  parameters,  and  then  moves  in  steps:  one  step  up,  one  step  right,  two  separate 
steps  down,  eva  separate  dteps  left,  three  steps  up,  and  so  on.  At  each  step,  the 
correlation  coefficient  is  calculated  as  described  in  section  3.2.  CC  is  thus  calculated 
at  each  position  along  the  spiral.  After  each  calculation  it  is  « ■■ninail  to  see  if  it 
is  any  larger  than  any  previous  CC  value  so  far  met  during  this  search,  end  if  so,  the 
new  CC  is  stored,  together  with  the  colt^n  and  row  ambers.  Dp  to  nine  previous 
values  in  those  arrays  arc  retained,  and  provided  CC  is  greater  then  0.3,  the  cumber 
of  the  current  ring  is  recorded,  at  all  times  retaining  the  number  of  the  ring  containing 
the  highest  value  of  CC  so  far  fo*od. 

Section  3.3  described  the  CC  surface,  which  in  pr»ct.ice  consists  of  a  hill, 
whose  peak  is  the  position  of  largest  CC  value,  rising  from  s  rough  plain.  The 
spiral  search,  starting  from  som  place  not  too  distant  from  the  kill,  cea  then  he 
imagined.  As  the  spiral  search  proceeds,  each  turn  rises  higher  up  the  flank  of  the 
bill,  until  eventually  the  peak  is  passed.  Thus  if  a  reasonably  high  CC  value  is 
met  and  then  no  higher  value  is  met  on  the  next  two  rings,  it  is  s  reasonable  wtuption 
that  the  peak  has  been  passed.  If  however,  no  CC  value  as  high  ss  say  0.5  has  keen 
met,  it  is  uncertain  that  the  peak  has  been  found,  and  the  search  is  allowed  to  centime 
until  all  the  rings  have  been  collated. 

As  described,  the  model  is  open  ro  objections  -  i:  is  possible  to  postulate  hills 
with  double  peaks,  or  oces  which  are  surrounded  by  s  high  plain  level.  It  ia  the 
essence  of  a  suitable  gep  that  it  does  not  have  such  behaviour.  Ideally  even  at 
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different  rises  of  yeu  or  with  different  son  elevation  etc.  the  CC  surface  should  show 
a  single,  fairly  steep  and  tall  hill  rising  froa  a  low  plain.  Cepe  Aick  do  not  shoo 
such  desirable  behaviour  should,  if  possible  he  manned  froa  the  list,  as  Che  fact 
becoaes  evident. 

In  practice  it  is  son* tines  foimd  that  a  gep  is  correctly  relocated  with  a  pash  CC 
value  less  than  0.5  and  in  such  cases  so as  confutation  tiae  rill  have  been  wasted  by 
continuing  the  search  unnecessarily,  however  it  is  not  possible  to  use  the  value  of  CC 
as  an  indication  of  true  relocation,  ae  experience  shone  that  values  froa  about  0.3 
upwards  can  correspond  to  correct  location,  and  peaks  of  up  to  0.6  or  worn  are  aenetiaea 
incorrect,  due  to  a  chatce  correspondence . 

Experience  has  shown  that  the  search  rules  described  above  uorfc  very  successfully, 
providing  that  a  little  cam  has  been  taken  in  the  choice  of  geps.  however,  gepe  which 
am  nis located  by  10CGCP  tend  to  be  those  difficult  to  locate  naaeally  ae  first 
identification,  so  in  practice  them  is  a  natural  bias  towards  'good*  geps.  The  aajor 
source  of  incorrect  locations  is  when,  for  sane  reason  such  s  a  change  of  ground  use  or 
the  gep  has  been  veiled  by  a  light  haze,  no  location  is  possible,  la  such  cases  IflCBg 
my  detect  a  position  where  CC  is  for  chance  reasons  greater  than  0.3,  and  this  madras 
location  is  returned  and  accepted.  Such  aislocations  can  usually  he  detected  ad 
eliminated  by  the  main  program  later  an. 

When  the  search  has  concluded,  some  results  are  displayed  on  the  user's  terminal. 
These  include  the  location  itself,  as  column  and  row,  together  with  the  peak  CC  valwe, 
hut  include  in  addition  a  table  of  CC  values  for  the  5*5  square  of  iaage  locatinan 
surrwding  that  peak,  and  the  (19  to)  ten  next  highest  CC  values  found  during  the 
search.  These  results  enable  the  user  to  visualize  the  CC  surface  and  the  approach 
of  the  spiral  to  the  correct  location.  Bich  of  this  infomatina  has  i sen  emitted  from 
a  recent  version  of  The  program,  as  sane  users  am  not  interested  in  it. 

7.1.)  Location  to  fractious  of  a  pixel 

Inspection  of  CC  surfaces  as  described  atom  shewed  that  of  tan  a  value  adjacent 
to  the  peak  was  only  slightly  lass  high,  indicating  that  the  'true*  location  was 
probably  between  the  two  integer  locations.  This  suggested  the  possibility  that  a  mam 
.curate  value  could  be  determined.  The  art  hod  chosen  was  to  start  from  the  integer 
location  of  maximm  CC  and  to  spiral  outwards  in  riags  of  0.1  pimel  rathe-  than  the 
usual  1.0  pixel.  To  perform  the  calculations  of  CC,  a  grid  of  pixel  values  at  intervals 
of  0.1  columi  and  row  are  needed,  aad  these  am  obtained  by  linear  interpolation  from  the 
adjacent  two,  or  more  usually,  four,  pixel  values.  This  method  enables  the  *  true* 
location  to  be  obtained  to  the  nearest  0.1  pixel.  Linear  interpolation  is  chosen  as 
probably  the  quickest  nr t bod.  and  although  clearly  sot  mathematically  the  nut  accurate 
is  probably  acre  than  sufficient  for  the  purpose,  consider  Lag  the  me  turn  of  the  overall 
system.  0.1  pixel  corresponds  to  about  6  u  in  the  scan  directum  and  S  a  in  the  travel 
direction. 

L00GCP  is  instructed  whether  to  search  to  an  integer  or  a  fractional  location  by 
as  ana  of  an  action  indicator  parameter.  A  value  of  2.1*  indicates  an  integer  search 
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«h«*M  a  value  of  1.0  (actually  any  value  not  within  the  range  of  about  1.5  to  2.5) 
indicates  that  a  fractional  search  is  required,  the  calculation  of  values  to  0.1  pixel 
spacing  is  in  principle  similar  to  that  that  for  integer  pixels.  The  creation  of  the 
spiral  is  done  by  a  series  of  steps,  as  before,  and,  as  for  the  integer  case,  is  stopped 
when  the  peak  CC  value  has  been  passed. 

7.2  Subroutine  MOE 

Since  tend  sat  satellite  orbits  are  neither  polar  nor  equatorial,  the  path  of  tike 
satellite,  or  rather  its  ground  nadir,  trill  make  some  non-zero  angle  with  true  north, 
the  transformations  described  in  thi  Beport  are  made  to  British  Kational  Grid,  the 
northings  of  uhich  are  only  true  uorth  at  cue  longitude  (2  degrees  west)  „  For  any  innge. 
the  centre  of  which  is  specified  by  a  path  and  row  uunber,  there  will  be  a  difference 
angle  between  satellite  path  and  grid  northing.  It  is  Che  purpose  of  sdbroutL  e  TWQF 
to  calculate  this  difference  angle  for  any  specified  path  and  raw  ssder. 

Inferring  to  the  net  bod  of  operation  of  the  overall  program,  once  the  first  gep 
has  been  located  attempts  are  wade  to  locate  other  geps  at  some  distance  f ran  die  first. 
The  approximate  locations  of  these  are  calculated  from  (amongst  other  things)  knowledge 
of  the  difference  angle  between  the  satellite  path,  (w  the  y  direction  of  the  usage) 
and  the  nap  northing.  nte  search  for  geps  is  limited  by  subroutine  LOCGCP  to  41  rings, 
re  up  to  41  pixels  in  all  directions.  Since  the  location  nay  be  iq  lo  say  1 00  km 
distant,  which  is  a  maximal  of  1 00  *  16  -  1600  pixels  away,  it  is  necessary  that  the 

difference  angle  be  known  to  an  accuracy  of  better  than  41  in  1600  or  about  1.4  degrees. 

It  is  therefore  important  that  1HCIF  provides  a  result  which  is  less  than  1.4  degrees 
in  error,  and  to  do  this  for  the  range  of  path  and  tow  numbers  applicable  to  all  images 
of  Britain. 

Is  noted,  the  difference  angle  consists  of  two  elements,  one  being  the  angle 
between  the  satellite  path  and  true  north,  and  the  other  the  difference  between  true 
north  and  Grid  north,  loth  of  these  angles  vary  with  latitude  and  in  dm  case  of  the 
second  with  Icegitude  also.  Since  the  insge  location  is  provided  as  a  path  and  row 
number,  it  is  necessary  firct  to  consort  these  to  longitude  and  latitude. 

The  operation  of  TMQF  is  therefore  the  following  angnence: 

(i)  The  latitude  and  longitude  of  the  centre  of  the  image  are  calculated  from 

the  Landsat  path  and  row  landii  r , 

(ii)  the  satellite  heading,  exprassed  as  the  angle  between  its  path  and  true 

north  is  calculated  for  that  latitude, 

(iii}  the  difference  angle  between  true  and  Kational  Grid  north  is  calculated 

for  that  latitude  and  longitude  (this  angle  is  referred  to  as  the  convergence) , 

(iv)  the  rotation  angle  is  calculated  by  addition  of  (ii)  aod  (iii). 

Items  (i)  to  (iii)  will  now  be  asidertd  in  nore  detail. 
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7.2.1  Ulculitioa  of  Utifte  loatiffc  fron  path  gj  row 

The  liiMt  ’Worldwide  hfcmct  Systeu*  of  paths  awl  row  is  described  is  Onf  2. 
■e  levant  details  iadak: 


(a>  There  are,  in  total,  251  evenly  spaced  paths,  iltkn^k  these  are  not 
traversed  in  newpence  by  the  satellite,  nevertheless  it  is  convenient  to  anAer  than  in 
a  sequential  nanaer.  The  separation  of  paths  is  therefore  340/251  -  1.4343  degrees, 
each  path  being  identical  to  its  ufighhnar  bet  displaced  in  longitude  by  this  aennt 
westward  fron  its  predecessor.  The  location  of  the  paths  is  not  defined  exactly  in  Kef  2 
in  terns  of  longitude,  although  it  is  noted  that  path  1  is  “the  first  track  that  includes 
mainland  Korth  4nerica“.  however,  candles  of  the  relationship  between  path/ roe  and 
longitude /latitude  are  provided  and  free  these  it  is  possible  to  calculate  the  longitude 
at  which  any  path  crosses  the  equator.  In  particular,  path  251  (which  is  the  sane  as 
path  'zero* )  is  found  to  cross  the  equator  at  44.05  da  gras  a  west.  (This  and  nil  other 
longitudes  in  this  Report  are  referred  to  the  Greenwich  neridian.) 


(b)  Each  conplete  orbit  takes  102  ui  water  (aore  exactly  4134  seconds)  and  iuegp a 
are  centred  on  locations  25  seconds  apart.  Thus  there  would  be  248  complete  iaages  in 
one  complete  orbit  (an re  exactly  247.84)  although  ouly  119  of  these  row  locations  are 
used.  The  119  are  recorded  os  the  aorth-south  portion  of  the  orhit  uhich  is  the  snalit 
side  of  the  earth,  sod  path  saber  40  is  at  the  equator,  how  1  is  the  next  northerly 
row  used,  being  at  80  degrees  1 .2  wiuures  north,  and  rev  1 19  the  no*':  southerly,  at  tha 
naan  latitude  south,  low  59  is  centred  on  a  location  25  seconds  before  the  equator  is 
reached,  and  so  on. 

Fran  this  intonation  it  is  possible  to  calculate  geographical  location  fron  ’"th 
and  row  ninhrrs. 


(i)  Latitude 

If  (Ik  orbit  of  the  satellite  is  as  a  on  i  to  be  a  circle,  it  nay  be  divided  into 
248  equal  sections,  corresponding  to  adjacent  row  angers.  Since  row  40  is  at  the 
equator,  any  row  BOH  i*  located  on  the  orbit  at  a  fractioa  (40-80H)/42  of  o  quadrant  fron 
die  equator,  sc 

latitude  of  B0U  on  the  orhit  •  p  *  *40-8010/124  radians. 

This  is  the  latitude  on  the  orhit.  which  is  at  an  iacliaatioa  aagle  A  to  tha 
equator.  Uben  the  orhit  is  related  to  the  ear*Vs  surface  tha  earth  latitude  Lat  is 
glaeo  by  the  expression 

sia(Lat)  -  sir.(A)  .  nia(pi  -  (4O-80U)/124) 

this  calcclatioo  being  for  a  sphere  rather  tkso  a  wore  caqplex  shape. 

|  V 

For  ao  ellipsoid 

tan(Lal_/  «•  e  *  a  »  tao(Lrtd) 
where  Late  is  the  geocentric  latitude  (Lat  above) 

Latd  is  the  geodetic  latitude  and 

is  ti  e  rtio  of  die  aiaor  to  the  naj«  r  axis 
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to  cross  the  motor  at  64.05 


Each  accewia  path  ia  di  ay  laced,  in  lagindc  only,  by  340/251  degrees  aaataai 
For  mj  giaca  path  therefore,  the  path  u^tr  am  be  Ml  tip  lied  by  340/251  to  provide 
the  eastward  ioacitede  shift  for  that  path,  relative  U  path  sero. 

Oar'icriat  aey  oaa  path,  the  loacitade  at  roe  aaher  BOB  ia  aa  socle  La  farther 
east  thee  at  row  aaftrr  40  (the  ages  tor)  where 

taaCLsa)  *  cos  (A)  .  taa(pi  *  (40-000/124) 
aad  this  expression  is  osed  to  evaluate  the  relative  laginde  for  aey  gives  roe  aA 

These  longitude  calculations  base  as  rased  that  the  earth  is  stationary,  whereas 
fact  it  is  rot  at  irg  is  relation  to  the  sea,  and  brace  the  satellite's  sai  syaihroanwa 
orbit,  at  the  rate  of  one  revelation  ia  24  boors.  For  each  increa.'e  is  row  wader  as 
interval  of  25  eemnds  baa  elapsed  dsriag  which  tier  the  earth  lies  rotated  by 
340  *  25/(24  v  3400)  decries,  U  0.10414  decree  per  raw.  It  is  therefore  accessary 
adjust  the  calculated  loacitade  by  this  amt  for  each  row  evader  differ  vary  free  th 
equatorial  row,  washer  40. 


it  for  each  row 


sax  ary  to 
fros  the 


steps  MceMay  for  the  calcelatiee  of  the  loacitade  occur 
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There  is  seen  to  be  fair 
be  doe  to  Such  c eases  as  si 


the  rill 


r as ideal  differences 


any  he  doe  to  each  causes  as  using  a  different  valse  for  e  (several  slightly  different 
'raises  are  in  one),  a  different  value  for  the  satellite  orbit  inclination,  or  to  sae  of 
a  suplificd  (spherical)  earth  nodal  for  soar  of  the  expressions. 

The  iaecr  supplied  for  any  particular  path  and  row  ia  nsnally  displaced  fna  its 
aoaiael  position  (sew  Kef  2) ,  often  by  e  larger  distort  then  the  above  err  are  ia  the 
calculated  locatiau.  Fx—netion  of  fifteen  i—y  shows  that  the  centre  am  be 
displaced  f roa  the  anaiwal  position  by  op  to  0.47  degree  in  Tnaagjtwde  aad  0.21  degree 
ia  latitude,  these  figures  including  tbs  shove  calculation  errors  of  abou 
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and  0.02  degree  respectively.  Tints  tie  calculated  longitude  ad  la*dtudt  are  adequately 
accurate  for  their  purpose.  The  effect  of  wee  error  in  location  of  t&s  centre  of  the 
inage  is  considered  later.  (If  ar, count  is  taken  of  the  satellite  in  ana,  it  is  fond 
that  the  longitude  displeceaents  for  Landaet  1  are  ia  the  range  0.337  to  0.404  degree, 
and  for  Lead* at  2.  -0.179  to  0.230  degree.  It  thus  eeean  that  Lands*  2  oca^ind  ei 
orbit  roughly  0.4  degree  eastward  of  Loudest  1  daring  the  relevant  .line  period.) 

7.2.2  Pete  mi  nation  of  the  satellite  heading 

The  satellite's  beading,  the  aqlc  he  mean  its  path  and  geodetic  north  at  a  given 
latitude  9  is  determined  with  the  aid  of  the  expression 

cos  A  •  cos  I  .  cos  C 

where  A  is  the  inclination  of  the  satellite  orbit  to  the  equator, 

B  is  the  latitude, 

C  is  rhc  inclination  of  the  satellite  orbit  to  that  latitude  parallel. 

The  inclinations  of  the  orbits  of  the  ladsat  series  are  provided  in  lef  2  ud 
measurements  of  the  orbits  are  recorded  in  lef  II.  Since  the  Lamdnnt  missions  are 
intended  to  provide  e  consistent  operational  system,  the  orbit  inclinations  do  not 
wary  greatly  mad  a  figure  of  99.1  degrees  may  he  used  for  all  of  the  satellites.  Since 
the  direction  of  satellite  travel  is  not  relevant  here,  the  angle  in  the  first  quadrant, 
80.9  degrees,  may  be  used  for  convenience.  Due  '„o  the  consistency  of  all  Landsat  orbits 
this  angle  will  in  general  be  in  error  by  appreciably  lesa  than  0.1  degree  and 
consequent  errors  in  C  should  he  less  than  0.1  degree  for  the  latitude  of  Britain. 

After  application  of  the  above  equation  at  the  appropriate  latitude,  the  angle  C 
is  converted  to  an  inclination  with  respect  to  the  local  meridian  rather  than  the  local 
parallel. 

7.2.3  Difference  between  true  north  and  Grid  north,  or  cemuergence 

The  British  Satiooal  Grid  refereue*  system  used  for  this  nosh  is  an  orthogonal 
system  superimposed  oa  a  transverse  Mercator  projection.  The  central  mtridin  of  the 
Grid  (2  degrees  west)  is  the  only  northing  line  which  ia  true  north,  all  other  northings 
having  some  convergence  from  true  north.  This  angle  any,  according  to  lef  12,  he 
approximated  by  the  expression 

convergence  *  (Lon  ♦  2)  .  sin(Lat) 

where  Lon  cad  Lat  are  the  longitude  and  lstitmle  ia  degrees.  This  expression  was 
coof  trued  by  comparison  with  a  table  of  values  calculsted  for  this  perpoas  and  is  in 
error  by  less  than  0.004  degree  for  ell  areas  which  are  normally  mapped  oa  the  British 
hatioaal  Grid. 

Further  increase  in  accuracy  is  not  justified,  es  the  present  expression 
contributes  negligible  error  to  the  subroutine. 


7.2.4  Accuracy  of  w>«witii»  TIMIf 

Ikt  owrall  acovicy  of  oftntin  of  auhroutine  Tlflf  is  4»t«raiw 4  by  tkt 
accuracies  of  its  several  cosstitoeet  ports.  These  have  hsss  discusaod  is  the 
appropriate  sebeectioss  and  css  nos  he  rn^iaii. 

Considering  first  the  longitude  ad  latitude,  a  stWy  of  f  if  toes  issues  detailed 
is  Table  1  shows  that  the  ceatre  of  as  image  asp  he  as  asth  ss  0.47  degree  ad 
0.21  degree  •'ey  froa  its  calculated  positios.  Swch  dif foresees  sap  he  shsea  to  give 
rise  to  errors  of  up  to  0.3S  degree  and  0.13  degree  is  the  result isg  differesce  angle. 
(These  figures  are  only  approximate  as  the;  ear;  with  latitude.) 

Io  addition  there  are  possible  errors  is  satellite  headiag  of  up  to  0.1  degree 
and  is  convergence  of  0.004  degree. 

Thus  when  all  errors  are  takes  into  account  it  appears  that  the  subroutine  IdTIF 
should  be  sble  to  provide  angles  shidk  are  is  error  by  less  than  abcut  0.5  degree. 

A  coup  arisen  has  been  ssde  between  the  rotation  aglw  provided  by  IMQ1  and  the 
actual  rotation  angles  de  term  sod  from  the  relevset  treeef onset  ins  matrices.  Fifteen 
images  of  widely  distributed  portions  of  the  CK  were  examined.  Each  image  provides  two 
rotation  angles,  one  of  which,  the  rotation  of  the  y-azis  of  the  image ,  is  due  to  the 
causes  described  in  this  section.  The  other  angle,  the  rotation  of  the  x-axis  of  the 
inapt,  contains  additional  smell  effects  doe  to  asy  yaw  amgle  of  the  satellite. 


Table  ! 

Differences  between  calonlatod  aad  actual  values  for  images  (degrees) 


Fsth 

- 

Centre  disy 
Longitude 

latitude 

yaagle 

Calculated 

Difference 

214 

23 

0.472 

0.013 

11.748 

11.891 

0.143 

217 

24 

9.041 

0.209 

13.372 

13.150 

-0.222 

218 

23 

01441 

-O.OiS 

13.945 

14.180 

0.215 

218 

24 

0.403 

0.041 

14.234 

14.273 

0.037 

219 

24 

0.141 

-0.042 

15.284 

15.394 

0.110 

219 

25 

0.140 

0.035 

15.194 

15.455 

0.241 

220 

22 

0.239 

0.014 

14.140 

14.404 

0.244 

220 

23 

0.204 

0.015 

14.297 

14.449 

0.172 

220 

24 

0.339 

0.058 

14.441 

14.519 

0.058 

229 

25 

9.337 

0.048 

14.843 

14.555 

-0.288 

221 

18 

-0.042 

-0.080 

17.015 

17.273 

0.258 

222 

19 

0.144 

-0.144 

18.147 

18.589 

0.442 

222 

20 

0.444 

0.059 

18.337 

18.452 

0.315 

224 

19 

0.184 

-0.108 

20.901 

21.035 

0.134 

224 

20 

-0.179 

-0.018 

20.804 

21.040 

0.254 

Then  when  the  angles  calculated  by  euhtwurisa  UMCU  were  cospered  with  the 
oe  repaired  by  the  in  age  y—exis  the  everagt  difference  found  was 


rotatii 


0.21 
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tbr  largest  difference  0.44  degree.  TWu  figures  are  consistent  with  the  expected 
accuracy  of  the  submit in*. 

Summarising,  the  overall  accuracy  of  subroutine  UMI  appears  to  he  consistent 
with  the  calculated  accuracy  of  about  0.5  degree,  adage  ate  for  its  intended  use,  which 
lrrJi  a  l ini t  of  1.4  degrees. 

It  should  he  eaphasised  that  the  use  of  subroutine  tMOf  is  only  appropriate  to 
the  Sritish  Katioual  Grid  ays'an  and  for  regions  sores l ly  napped  to  that  Grid.  This 
does  not  include,  for  exnaple.  any  part  of  Ireland. 

The  mathenatical  expression*  used  are  in  general  not  rigorous  (for  exanple  sane 
earth  node  la  are  based  on  a  sphere  rather  then  an  ellipsoid),  but  are  seen  to  he 
adequate ly  accurate  representations. 

7.3  Subroutine  HSTSCh 

Subroutine  ISTStS  ha*  the  purpose  of  calculating  a  first-order  transformation 
matrix  from  a  set  of  pairs  of  locations  in  the  two  reference  frames:  image  and  Sational 
Grid  coordinate  system*.  This  subroutine  is  in  fact  adapted  Iran  an  existing  program 
NUIIX*  so  it  will  not  be  described  in  detail  here,  lusuir  it  now  contain*  some 
modifications,  as  follows. 

arSCI  has  been  altered  so  that  it  reads  the  corresponding  sat*  of  image  locations 
and  nap  coordinates  fran  appropriate  arrays  via  a  COMBS  block,  rather  then  by  obtaining 
then  from  a  data  file.  In  addition,  it  ones  a  small  amount  of  discrimination.  If  any 
gep  has  a  residual  (determined  by  a  previous  HalSCa  calculation)  of  greater  than  as 
wnount  apecified  by  one  of  the  subroutine  parameters,  or  has  a  correlation  coefficient 
of  less  than  zero,  that  gep  is  not  included  in  the  currant  calculation.  Initially  all 
residuals  are  set  to  zero  by  the  main  program,  so  the  first  N4XSEB  operation  uses  all 
avai table  geps.  A  correlation  coefficient  of  less  than  aero  iaplics  that  the  gep  has, 
for  one  reason  or  another,  not  been  relocated:  in  fact  say  CC  value  less  than  0.3  will 
he  provided  to  HATSCB  as  less  than  zero. 

If  less  than  four  suitable  gepe  are  provided,  MUSIS  is  not  able  to  fwetion,  so  it 
provides  an  appropriate  error  aessage,  and  then  stops.  This  eight  occur  if  the  new  iange 
wire  largely  obscured  by  cloud,  so  that  few  gep*  could  he  relocated. 

Uhea  MITSCI  has  performed  the  least  ngaates  fitting,  it  then  prepares  a  new 
transformation  matrix  which  is  pass  ad  to  the  neia  program  through  the  CUHBS  block  Mlllll. 

During  the  course  of  its  operation,  MUSIS  provides  sons  output  to  th*  near 
terminal .  concerning  both  the  details  of  the  individual  gep*,  their  locations  and 
residuals,  and  also  details  of  the  now  transformation  matrix,  including  the  inplied 
pixel  sizes  and  the  relative  orientation  both  slang  and  perpendicular  to  the  satellite 
direction  of  travel.  This  data  output  has  been  omitted  from  a  recent  version  of  the 
program. 

MATSU  therefore  is  a  fairly  wopbi sticated  subroutine »  providing  only  a  first- 
order  trwsfornatiaa  matrix.  It  is,  however,  quite  adequate  for  in  GCf.FDt,  and 


■on  complication  is  sot  justified.  One*  a  sot  of  gepe  has  bene  relocated  using  CCF.FDD 
the  user  nap,  if  ha  wishes,  usa  a  program  CCfflT  to  obtain  a  first-  or  higher-order 
tranefvrmetiom  matrix,  selecting  such  gep  locations  as  ha  mgnirea. 

7.4  library  end  spates  suhrontinas 

7.4.1  Snhroutines  in  library  HDKJB 

Program  OCF.FDB  makes  see  of  several  subroutines  hold  in  a  library  mimed  HH.I1. 
The  HH.D  routines  used  are  listed  hare,  with  a  brief  description  of  their  action. 

CMH:  ton.  a  chip  filename  from  a  gep  identity. 

CLSOt:  closes  a  temporary  file. 

OIEX.  opens  an  image— format  file  for  reading,  using  the  filenane  aa  a  par  neater. 
OMPgt:  opens  a  temporary  file  fer  reading  or  writing. 
ailUL:  reads  the  toil  (if  any)  of  an  image-format  fila. 

TOCS:  outputs  to  the  specified  channel  (l  *  terminal)  details  of  the  data  and  certain 

tines.  This  allows  estimates  to  he  nede  of  the  it  apt  of  coagutor  processing  and 
input/output  tine. 

7.4.2  Subroutines  in  library  1HLIB 

One  is  made  of  none  subroutine*  held  in  the  inage-fornat  library  1WJ1.  These  on 
CLOSE I :  donee  «  image— format  file. 

miTI:  imitializes  data  in  «  CHIOS  block  used  by  OLD  routines. 

OIK:  opens  an  image— format  file  for  reading. 

7.4.3  Systen  subroutines 

The  program  GCP.FIHD  has  been  developed  for  use  on  e  Prime  730  cenpnter  and  is 
written  in  friar  F0BT1AX  IT.  It  makes  use  of  a  mtsdter  of  Prime  system  embrnu rinse, 
most  of  which  ere  of  a  familiar  nature,  such  es  Kid  or  tHIK.  lose  familiar  routines 
include: 


QjOSU:  closes  e  file. 

MDU:  'rewinds*  a  file,  v<  returns  to  the  start. 

TSKOfk:  presents  e  question  at  the  user  terminal  amd  sets  the  logical  value  HI  ee 

FALSE  according  to  the  response  YES  or  HO. 

It  ie  necessary  to  include  the  iaoert  STSOMMlms  in  the  declaration  section  of 
ell  progr—  md  subroutines  which  use  ouch  systen  subroutines,  to  'explain*  to  the 
caepiler  the  value  of  words,  used  as  psrmnters  ie  these  smbrou  times,  which  is  elude  the 
*«*  character. 

8  gESOUS  MED  COPCLCSIOB 

The  'product* of  GCP.FIHD  is  a  list  of  all  the  geps  uhich  hams  been  relocated 
within  e  given  image.  Ideally,  ell  gepe  which  lie  within  the  imp  will  he  relocated. 


33 


ad  to  *  high  accmracy.  In  practice,  Am  to  aaall  aulocn  ions  of  the  original  yya,  a 
accuracy  of  about  £0.4  pixel  it  to  ha  npoctad,  rather  tbaa  tba  0.1  pial  aypiraailj 
possible. 

Tba  alternative  to  wsiag  OCf.IUO  la  to  locata  aa  at?  aa  pouiklt  gtyt  aaully, 
which  ia  a  tfat  cowsuwing  labour  iatamim  process.  Tba  aoccaaa  or  othexwiae  nf  Ctf-MM 
ia  tharofora  to  ha  aaaaurad  by  where  ita  pfffonact  lioa.  whether  'perfect* ,  or  'butter* 
or  'worse*  tbaa  tba  aaaoal  athod. 

To  pnvidt  a  gwalitatiue  answer  to  tkia  foatioo,  it  aoold  bo  oacaaaaty  to  mold 
systaustically  tba  raaolta  of  aaa  of  Off.FDB  a  a  *ida  iafi  of  iat^t,  ad  no  rarb 
ayataaatic  r» cording  baa  k«  dooa.  Hjwaatr  tba  progran  baa  baaa  oaad  oa  oaay  diffcrta 
iasges  (over  100)  by  'b»  author  ad  aemral  otbar  uaers.  and  tba  follwiig  onaut 
nprtaat  tba  experience  fouad. 

CCr.FUD  ia  aery  acb  acre  cooaaoiaat  to  aaa  tbaa  tba  aaaoal  aatbodl  aid  ao  war 
would  now  cxmteaplate  aacoal  ralocatioa.  except  possibly  for  aa  iap  ahich  aaa  ao 
cloud  coaa rad  that  oaly  a  wry  aaall  amber  of  widaly  scattarad  gepe  could  ba  discerned. 
Tba  pragma  caa  cactaialy  ba  said  to  bo  vary  aacb  batter  than  aaaaal  aathoda. 

Tba  proportion  of  gepe  relocated  within  aa  inage  varies,  depending  oa  tack  things  as 
tba  season  of  the  year  cij^aiad  with  nfaan  tba  gep  chips  were  obtaiaad.  fcr  siailar 
re  in  oat,  the  proportion  relocated  approach!  ■  I0QX,  whilst  in  the  norst  caaas  attempted 

ahes;  ££Z  si  the  gep*  were  fsesd,  abas  tba  issge  was  abet  six  aeetbs  differest  in 

season  from  the  original  chips,  and  sane  snow  present  ia  sane  portions  of  tbs  scans. 

Tba  accuracy  of  ralocatioa  any  ba  judged  fron  tba  standard  deviations  of  tba  gepe. 

In  tba  best  cases  the  standard  deviation  nay  ha  lass  than  30  n  (about  0.4  to  0.9  piael) 
with  one  or  two  out  of  perhaps  ISO  being  bad  'outliers* .  la  the  uorst  situation,  the 
standard  deviation  nay  he  120  n  with  up  to  1 5-201  badly  relocated,  aaialy  because  they 
cannot  he  located  at  all. 

It  should  be  pointed  oat  that  heavy  cloud  cover,  or  * 'marines  in  Britain  poor 
atamsphsric  clarity,  reduce  the  value  of  inages,  ao  thane  results  refer  to  inagse  which 
appear  usable  after  a  wry  brief  visual  inspection. 

In  wary,  CCP.FDKD  provides  s  rapid  nrthod  of  relocating  gepe,  with  little  near 
effort,  nd  enables  large  guaatitias  of  jaggery  to  ha  process  ad,  nhich  would  not  kare 
Wan  contenplsted  with  nar  ml  net  hods.  Provided  that  the  inage  quality  ia  good,  tW 
process  has  a  high  success  rata  in  accurately  relocating  geps. 
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